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ABSTRACT 


An  analytical  approach  toward  numerical  calculation  of  the  three- 
dimensional  turbulent  boundary  layer  on  a  sharp  cone  at  incidence  under 
supersonic  and  hypersonic  flow  conditions  is  presented.  The  theoretical 
model  is  based  on  implicit  finite -difference  integration  of  the  governing 
three-dimensional  turbulent  boundary-layer  equations  in  conjunction  with 
a  three-dimensional  scalar  eddy- viscosity  model  of  turbulence.  Com¬ 
parison  of  the  present  theory  with  detailed  experimental  measurements 
of  the  three-dimensional  turbulent  boundary-layer  structure  (velocity 
and  temperature  profiles),  as  well  as  surface  streamline  direction 
(obtained  via  an  oil-flow  technique)  and  surface  heat-transfer  rate,  re¬ 
veals  good  agreement.  Effects  of  wall  temperature  on  the  three- 
dimensional  turbulent  bound  ary -layer  profiles  under  hypersonic  condi¬ 
tions  are  considered  relative  to  interpretation  of  hot  wall,  hypersonic 
wind  tunnel,  force  testing  as  it  relates  to  cold  wall,  atmospheric  re-entry. 
The  calculated  surface  upwash  angle  is  found  to  be  fairly  sensitive  to 
wall  temperature  effects  with  the  larger  values  of  the  angle  occurring 
with' the  hotter  wall. 
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SECTION  I 
INTRODUCTION 


Only  within  the  last  decade,  with  the  advent  of  large,  high-speed 
digital  computers,  has  the  problem  of  numerical  integration  of  the  two- 
dimensional  compressible  turbulent  boundary -layer  equations  using  a 
two-layer  (inner-outer)  eddy-viscosity  model  to  describe  turbulence 
become  feasible.  The  works  of  Smith  and  Cebeci  (Ref.  1),  Patankar 
and' Spalding  (Ref.  2),  and  Herring  and  Mellor  (Ref.  3),  as  well  as  the 
Langley  Compressible  Turbulent  Boundary  Layer  Symposium  (Ref.  4), 
are  representative  examples  of  the  current  state-of-the-art  in  numeri¬ 
cal  calculation  techniques  for  compressible  two-dimensional  turbulent 
boundary -layer  flows.  In  general,  the  two-layer  (inner-outer)  eddy- 
viscosity  law,  coupled  with  the  so-called  van  Driest  damping  for  the 
near-wall  region,  appears  to  be  entirely  satisfactory  and  accurate  for 
two-dimensional  turbulent  boundary-layer  flows  in  the  subsonic  and 
supersonic  regime.  More  recent  works  (Refs.  5-9)  indicate  that  the 
same  two-layer  eddy -viscosity  law  is  applicable  even  in  the  hypersonic 
regime. 

With  regard  to  the  three-dimensional  compressible  turbulent 
boundary -layer  problem,  little  work  has  been  done  to  present.  The 
three-dimensional  compressible  turbulent  boundary -layer  equations 
have  been  derived  by  Braun  (Ref.  10)  and  Vaglio-Laurin  (Ref.  11). 
Several  attempts  to  solve  these  equations  using  integral  techniques  in 
conjunction  with  the  so-called  small  crossflow  assumption  (which  un¬ 
couples  the  streamwise  momentum  equation  from  the  crossflow  mo¬ 
mentum  equation)  have  been  made  by  Cooke  (Ref.  12),  Smith  (Ref.  13), 
and  Bradley  (Ref.  14),  as  well  as  Braun  and  Vaglio-Laurin  cited  pre¬ 
viously.  These  integral  approaches  were  not  entirely  satisfactory  be¬ 
cause  of  the  difficulties  in  adequately  representing  the  crossflow  velocity 
profile.  However,  the  recent  work  of  Shanebrook  and  Hatch  (Ref.  15) 
employing  hodograph  models  for  the  crossflow  velocity  component  ap¬ 
pears  promising  toward  removal  of  the  above-mentioned  difficulty.  It 
should  be  mentioned  in  this  connection  that  Zakkay  and  Calarese 
(Ref.  16)  have  successfully  applied  an  integral  analysis  for  the  com¬ 
pressible  turbulent  boundary  layer  undergoing  both  adverse  pressure 
gradient  and  crossflow  along  a  plane  of  symmetry. 

The  paper  by  Hunt,  Bushnell,  and  Beckwith  (Ref.  17)  is  the  first, 
to  the  author's  knowledge,  to  apply  marching  finite -difference  .integra¬ 
tion  to  the  three-dimensional  compressible  turbulent  boundary-layer 
equations  for  swept  infinite  cylinders  under  hypersonic  conditions.  In 
this  work  three  different  approaches  to  the  formulation  of  eddy -viscosity 
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models  for  three-dimensional  turbulent  boundary  layers  were  con¬ 
sidered  with  two  of  the  formulations  being  used  in  the  numerical  solu¬ 
tions  and  one  being  completely  rejected  (the  independence  principle 
approach).  Of  the  two  acceptable  eddy-viscosity  models,  the  most 
plausible  formalism  (in  the  opinion  of  the  present  author)  is  the  in¬ 
variant  turbulence  approach  which  is  based  on  the  concept  that  the 
three-dimensional  eddy  viscosity  should  depend  only  on  the  properties 
of  the  turbulence  and  a  local  eddy  scale.  The  application  of  this  con¬ 
cept  to  a  mixing-length  model  suggests  that  the  eddy  viscosity  should 
be  a  scalar  function  independent  of  the  coordinate  direction.  Implicit 
in  this  approach  is  the  requirement  that  the  directions  of  the  three- 
dimensional  turbulent  shear  stress  components  are  the  same  as  the 
directions  of  the  corresponding  mean  velocity  gradients.  As  shown  by 
Bradshaw  (Ref.  18)  through  a  three-dimensional  kinetic -energy- of  - 
turbulence  approach,  the  directions  of  the  shear  stress  components 
are  not,  in  general,  exactly  the  same  as  the  directions  of  the  cor¬ 
responding  mean  velocity  components.  However,  in  the  near -wall 
region,  Bradshaw's  shear  stress  equations  reduce  to  the  scalar  eddy 
viscosity -mixing  length  model  discussed  above. 

Recent  analyses  by  Nash  (Ref.  19)  and  Cooper  (Ref.  20)  concerning 
three-dimensional  turbulent  boundary  layers  in  an  incompressible  flow 
indicate  that  use  of  the  above -discussed  invariant  turbulence  concept 
will  lead  to  acceptable  results  when  coupled  to  an  accurate  integration 
technique  for  numerical  solution  of  the  three-dimensional  turbulent 
boundary- layer  equations.  Nash  uses  an  empirically  modified  kinetic- 
energy -of-turbulence  equation  in  conjunction  with  the  assumption  that 
the  turbulent  shear  stress  acts  in  the  direction  of  the  mean  rate  of 
strain  of  the  mean  motion  to  compute  the  magnitude  of  the  turbulent 
shear  stress.  The  work  by  Cooper  applies  the  scalar  eddy  viscosity¬ 
mixing  length  model  to  the  problem  of  three-dimensional  incompressible 
laminar  and  turbulent  boundary -layer  flow  associated  with  a  rotating 
disk  in  an  infinite  fluid  otherwise  at  rest.  Cooper  reports  excellent 
agreement  between  theory  and  experimental  data  for  radial  and  circum¬ 
ferential  velocity  profiles,  displacement  thickness,  shape  factor,  and 
drag  attributable  to  skin  friction  for  both  laminar  and  turbulent  three- 
dimensional  boundary  layers.  It  should  be  noted  here  that  Cooper 
applied  a  modified  version  of  the  Cebeci  and  Smith  (Ref.  21)  linearized 
implicit  finite -difference  method,  whereas  Nash  used  an  explicit  finite- 
difference  scheme  which  has  recently  been  updated  and  improved  (see 
Ref.  22  for  clarification).  The  method  of  Nash  has  also  recently  been 
modified  to  include  the  effects  of  centrifugal  and  Coriolis  forces  for 
application  to  calculation  of  the  incompressible  three-dimensional  tur¬ 
bulent  boundary  layer  on  a  helicopter  rotor  as  reported  in  Refs.  23 
and  24. 
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The  present  document  reports  on  an  analytical  investigation  and 
development  of  a  finite -difference  calculation  technique  for  the  analysis 
of  the  three-dimensional  compressible  turbulent  boundary  layer  on  a 
sharp  cone  at  angle  of  attack  in  a  supersonic  or  hypersonic  flow.  The 
sharp  cone  geometry  at  angle  of  attack  is  of  obvious  importance  to 
aerodynamic  and  propulsion  engine  applications.  Specifically,  the 
governing  three-dimensional  turbulent  boundary -layer  equations  for  a 
compressible  flow  are  simplified  by  the  conical  nature  of  the  sharp 
cone  flow  field  and  numerically  integrated  on  a  digital  computer  utilizing 
a  marching  implicit  finite -difference  technique.  Three-dimensional 
turbulence  is  accounted  for  using  the  three-dimensional  eddy -viscosity 
invariant  turbulence  approach  discussed  above  in  conjunction  with  an 
inner-outer  mixing-length  formalism  from  two-dimensional  flows 
carried  over  to  the  three-dimensional  flow  of  present  interest.  The  in- 
viscid  conical  flow  field  about  the  sharp  cone  at  incidence  is  deter¬ 
mined  using  a  prior  documented  digital  computer  code  which,  in  turn, 
furnishes  the  outer-edge  conditions  for  input  to  the  present  boundary- 
layer  analysis.  To  assess  the  accuracy  and  applicability  of  the  present 
theory,  comparisons  of  the  calculated  three-dimensional  turbulent 
boundary -layer  profiles  (velocities  and  temperature)  are  made  with 
detailed  experimental  flow-field  surveys  of  the  three-dimensional  tur¬ 
bulent  boundary-layer  structure  at  various  circumferential  locations 
around  a  sharp  cone  at  incidence  under  supersonic  flow  conditions. 
Comparisons  between  the  present  theory  and  experimental  measure¬ 
ments  of  surface  heat  transfer  under  hypersonic  flow  conditions  are 
also  given  for  a  sharp  cone  at  various  angles  of  incidence.  Detailed, 
hypersonic,  three-dimensional,  turbulent  boundary -layer  profiles 
around  a  sharp  cone  at  incidence  are  presented  (although  there  are  no 
experimental  profiles  for  comparison).  Effects  of  wall  temperature  on 
the  three-dimensional  turbulent  boundary  layer  under  hypersonic  condi¬ 
tions  are  examined  in  some  detail  relative  to  interpretation  of  hot  wall, 
hypersonic  wind  tunnel  force  testing  as  it  relates  to  cold  wall,  atmos¬ 
pheric  re-entry. 


SECTION  II 

ANALYTICAL  ANALYSIS 


The  present  analytical  investigation  employs  a  three-dimensional 
turbulent  boundary-layer  analysis  coupled  with  a  three-dimensional  in- 
viscid  conical  flow  analysis  for  a  sharp  cone  at  incidence  in  a  supersonic 
or  hypersonic  stream.  Development  of  the  boundary -layer  analysis  is 
presented  below;  the  inviscid  analysis  utilizes  a  documented  digital  com¬ 
puter  code  which  is  described  briefly. 
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2.1  GOVERNING  THREE-DIMENSIONAL  TURBULENT 

BOUNDARY-LAYER  EQUATIONS 

The  present  analysis  employs  the  three-dimensional  compressible 
turbulent  boundary-layer  equations  in  terms  of  time -averaged  mean  flow 
quantities  as  derived  by  Vaglio-Laurin  (Ref.  11).  The  coordinate  sys¬ 
tem  is  taken  to  consist  of  geodesics  and  geodesic  parallels,  following 
Moore  (Ref.  25).  In  this  coordinate  system  the  body  surface  is  defined 
by  y  =  0  and  a  point  in  space  is  defined  by  the  distances,  x,  y,  and 
r(x)0  where  the  length  r(x)$  depends  explicitly  on  the  distance  x  and 
where  r(x)  has  the  dimensions  of  length.  This  coordinate  system  is 
especially  useful  in  the  analysis  of  flow  about  bodies  for  which  a  co¬ 
ordinate  x  can  be  defined  such  that  body  cross  sections  are  similar  for 
various  values  of  x.  The  quantity  r(x)  then  gives  the  variation  of  scale 
for  these  cross  sections.  For  the  sharp  cone  geometry  under  present 
consideration,  the  geodesic  coordinates  are  taken  to  be  the  cone  gen¬ 
erators  and  the  geodesic  parallels  are  the  circles  swept  by  the  meridional 
angle.  The  corresponding  length  function  r(x)  is  the  local  radius  of  the 
body.  See  Fig.  1  for  clarification  of  the  sharp  cone  geometry,  nomen¬ 
clature,  and  coordinate  system.  The  time-averaged  mean  velocity  com¬ 
ponents  are  taken  to  be  u,  v,  and  w  in  the  directions  of  x,  y,  and  0, 
respectively.  The  governing  equations  of  motion  are,  following  Vaglio- 
Laurin  (Ref.  11): 


CONTINUITY 


r*  (^nr)  +  h  (?Vr)  +  h  (?w)  =  0 


STREAMWISE  (x)  MOMENTUM 


_ <?u  _,.(5u  .  pw  (Ju 

P«  —  +  p  V—  +  —  31 


dx 


dy  «■  d  <f> 


p(w) 2  <9  r 
r  <9x 


d  r  <?u  sr-r-Tl 
dy  t  dy  J 


CIRCUMFERENTIAL  ( <j> )  MOMENTUM 


_ dw  <9w  p  vv  dw  p  u  w  dr 

Pu  HZ  ~  PV3F  +  —  36  +  37 


-1  d\)  d  ["  dw  _  •  ■; 

=  7^  3^37  ~  P  V"J 


NORMAL  (y)  MOMENTUM 


JE  =  0 

dy 


(1) 


(2) 


(3) 


(4) 
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ENERGY 


where 


__5H  _..9H  pwffl  d  f  (dH 

puj^  +  px'w+  =  37[W  + 


^  $-  H 


V  = 


v  + 


/  / 
P_V 

9 


(5) 


(6) 


and  the  usual 'expressions  for  the  mean  and  fluctuating  parts  of  the  de¬ 
pendent  variables  are  used;  e.  g. , 


P  =  9  +  9 


(7) 


Implicit  in  Vaglio-Laurin's  derivation  of  the  above  equations  are  the 
following  stipulations : 

a.  The  rates  of  change  of  the  mean  flow  properties  in  the 
x-  and  ^-directions  [0(1)]  are  smaller  than  the  rates 
of  change  in  the  y-direction  [0(6”^)]  by  an  order  of 
magnitude. 

b.  Mean  squares  and  products  of  the  turbulent  fluctuations 
are  0(6);  that  is,  the  turbulent  level  is  small.  The 
terms  involving  mean  squares  of  the  velocity  fluctua¬ 
tions  are  taken  to  be  negligible,  which  is  valid  for  high 
Reynolds  number  flows  with  a  zero  or  favorable  pres¬ 
sure  gradient. 

c.  The  time  average  molecular  transports  are  approxi¬ 
mated  by  those  pertaining  to  the  mean  flow  properties; 
indeed,  even  the  latter  are  negligible,  except  very  near 
the  wall,  compared  with  terms  involving  the  turbulent 
transports. 

Also  implicit  in  the  above  equations  is  the  requirement  of  an  inviscid 
conical  flow  field  which  leads  to  the  term  9p/9x  =  0;  a  full  discussion  of 
the  inviscid  conical  flow  field  about  a  sharp  cone  at  incidence  in  a  super¬ 
sonic  or  hypersonic  flow  will  be  given  later  in  this  section. 


The  energy  equation  (5)  is  defined  in  terms  of  the  mean  stagnation 
enthalpy 


H  =  h  + 


u  2  +  w2 
2 


(8) 


If  the  two  momentum  equations  ( 2)  and  (3)  are  multiplied  by  u  and  w, 
respectively,  and  added,  one  obtains  an  equation  for  (u2  +  w2)/2  which 
is  simplified  by  the  fact  that  the  curvature  terms  (those  involving  9r/9x) 
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vanish.  If  then  this  equation  is  combined  with  the  energy  equation  (5), 
one  obtains  the  following  energy  equation  in  terms  of  the  mean  static 
enthalpy: 

ENERGY 


o 


This  is  the  form  of  the  energy  equation  used  in  the  present  analysis. 

If  subscript  w  denotes  wall  and  subscript  e  denotes  outer  edge  of 
the  boundary  layer,  the  associated  boundary  conditions  on  the  above 
defined  equations  are 


MOMENTUM 


ENERGY 


y  =  0  :  u  =  v  =  w=  u  v'  =  v'w'  =  p  'v '  =  0 

lim  y  -*  oo  :  u  -*  Ue,  w  -*  Wg 

u  'v '  0,  v  'w '  0 ,  p  'v '  -*  0 


(10) 


y  =  0  :  h  =  hw,  v'h'  =  0 

-  —  (11) 

which  reflect  the  requirements  of  no  slip  and  no  homogeneous  mass  in¬ 
jection  (suction  or  blowing)  at  the  wall  as  well  as  a  prescribed  constant 
wall  enthalpy.  The  normal  momentum  equation  (4)  reveals  that  the 
static -pressure  variation  across  the  boundary  layer  is  negligible,  and 
hence  the  static  pressure,  p(<t>),  is  regarded  as  an  external  input  to  the 
boundary -layer  analysis  from  a  separate  inviscid  analysis.  The  outer- 
edge  velocities,  Ue  and  We,  as  well  as  the  outer-edge  static  enthalpy, 
he,  must  be  determined  from  the  inviscid  analysis  consistent  with  the 
imposed  static-pressure  distribution. 

The  gas  model  adopted  for  the  present  study  is  thermally  and 
calorically  perfect  air  having  a  constant  specific  heat  ratio  y  =  1.40  and 
obeying  the  equation  of  state 

P  =  pRT  (12) 
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where  R  =  1716  ft^/sec^-°R.  Hence,  under  this  assumption, 

£  =  CpT  ( 13) 

where  Cp  =  6006  ft^/sec^-oR.  The  laminar  viscosity,  M,  is  taken  to 
obey  Sutherland's  law 

fl  _  ^ref  -r  198,6  /  T  \  (14) 

niel  ~  T  +  198.6  \Tre{/ 

where  T  must  have  units  of  °R  and  subscript  ref  denotes  a  reference 
condition.  The  laminar  Prandtl  number,  Pr,  is  taken  to  be  a  constant 
value  of  0.  71  across  the  entire  boundary  layer. 


2.2  TURBULENT  TRANSPORT  LAWS 


Before  Eqs.  (1),  (2),  (3),  and  (9)  can  be  solved,  expressions  must 
be  supplied  for  the  Reynolds  stress  or  turbulent  shear  terms  in  the 
momentum  equations  and  the  turbulent  flux  of  static  enthalpy  in  the 
energy  equation.  The  approach  used  in  the  present  analysis  is  to  model 
these  terms  as  functions  of  the  mean -flow  variables  following  the 
analysis  by  Hunt,  Bushnell,  and  Beckwith  (Ref.  17). 


The  concept  that  the  Reynolds  stress  in  turbulent  flow  is  propor¬ 
tional  to  a  momentum  exchange  coefficient  times  the  mean-flow  velocity 
gradient  normal  to  the  surface  is  well  known  and  commonly  used  in  tur¬ 
bulent  boundary -layer  analyses.  This  concept  is  based  on  an  assumed 
analogy  between  the  so-called  eddy  viscosity  and  the  molecular  viscosity. 
The  total  shear  components  in  the  streamwise  (x)  and  circumferential  ($) 
directions  are  written  as 


dll  —-T-7 

T  =  U  -r—  —  P  U  V 

*  p  dy  v 

du  du 

=  *Ty  +  '*Ty 

(15) 

dw  _  —7—7 

dw  dw 

(16) 

V  =  ^  -  p  vw 

^  dy  +  t<t>  dy 

where  the  eddy  viscosities  ex  and  in  the  x-  and  ^-directions,  respec¬ 
tively,  might  in  general  be  different.  Sioce  the  total  resultant  shear 
must  be  a  vector  quantity,  its  magnitude  is  written  as 


H  +  € 


W-mW  (i7> 
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In  a  similar  manner  the  total  heat  flux  in  the  energy  equation  (9)  be¬ 
comes 

k  <?h  _  ■  7  k  <Jh  k  3h 

^  C  dv  ^  V  C  <?v  +  C  <?v  (18) 

where  k  is  the  so-called  eddy  thermal  conductivity. 


2.3  EDDY-VISCOSITY  MODEL 


The  simplest  approach  to  the  formulation  of  models  for  the  Reynolds 
stress  is  based  on  Prandtl's  mixing-length  hypothesis.  For  two- 
dimensional  flow,  this  hypothesis  states  that 


'-ish 


(19) 


The  turbulent  shear  and  eddy  viscosity  are  then 


’  lurh 


-p  U  V 


Ou  o 

%■ 


dn 


from  whence 


f 


pf 2  & 
t  cry  I 


(20) 


(21) 


The  quantity  is  termed  the  mixing  length  and  is  some  characteristic 
length  related  to  the  size  or  scales  of  eddies  responsible  for  the  flux  of 
momentum  in  the  y-direction.  Although  the  details  of  such  a  transfer 
mechanism  are  not  well  understood,  the  basic  concept  gives  satisfactory 
results  even  at  hypersonic  Mach  numbers  in  the  presence  of  heat- 
transfer  and  pressure  gradients  as  mentioned  in  the  review  of  turbu¬ 
lent  boundary-layer  literature  given  in  Section  I. 


The  turbulent  boundary -layer  studies  by  Prandtl  (Ref.  26)  and 
Glushko  (Ref.  27)  are  based  on  the  concept  that  the  eddy  viscosity  should 
depend  only  on  the  properties  of  the  turbulence  and  a  local  eddy  scale. 
Extending  this  concept  to  the  three-dimensional  case  of  present  interest 
suggests  that  the  eddy  viscosity  should  be  a  scalar  function  independent 
of  the  coordinate  direction.  Accordingly,  the  components  of  the  Reynolds 
stress  are  written  as 


—  /  / 

—  p2  <?G  <?u 

pt •  575^ 

7turb,  x  s  ~P  u  v  = 

-  /  / 

-2  <*G(9w 

rlurb,  <}>  =  V  "  " 

pL  dyW 

(22) 

(23) 
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so  that  the  eddy  viscosity  is 


(24) 


To  determine  the  scalar  function  G,  Eq.  (24)  is  substituted  into  Eq.  (17) 
with  the  result  that  the  magnitude  of  the  total  resultant  shear  becomes 


r  a 


(25) 


Now,  by  analogy  with  Eqs.  (20)  and  (21)  where  the  velocity-gradient 
function  for  the  turbulent  shear  component  is  repeated  in  the  eddy- 
viscosity  expression,  the  relationship  for  the  determination  of  the 
G  function  becomes 


<9G 

57 


which  results  in  the  final  form  of  the  eddy  viscosity 


(26) 


(27) 


% 

Hunt,  Bushnell,  and  Beckwith  (Ref.  17)  call  the  above  representa¬ 
tion  of  the  Reynolds  stress  the  invariant  turbulence  model.  The  key 
assumption  in  this  model  is  that  the  eddy  viscosity  is  a  scalar  function 
independent  of  coordinate  direction.  In  a  recent  analysis  of  three- 
dimensional  incompressible  turbulent  boundary -layer  flows  using  the 
kinetic -energy-of-turbulence  approach,  Nash  (Ref.  19)  has  advanced 
arguments  that  the  turbulent  shear  stress  is  likely  to  act  in  the  mean 
rate  of  strain  direction,  defined  by  the  components  of  the  mean  velocity 
gradient  vector,  so  that  his  closure  equation  is  simply,  in  the  nomen¬ 
clature  of  the  present  analysis, 

rturb,  x  rturb,  tf> 

du  8w  (28) 

By  dy 

The  same  equation  results  from  the  scalar  eddy -viscosity  model  pre¬ 
sented  by  Eqs.  (22)  and  (23)  above.  On  the  other  hand,  Bradshaw 
(Ref.  18)  has  derived  a  set  of  differential  equations  for  the  two  com¬ 
ponents  of  the  turbulent  shear  stress  based  again  on  the  kinetic -energy- 
of-turbulence  approach  but  permitting  the  turbulent  shear  stress  vector 


9 


AEDC-TR-72-66 


to  deviate  from  the  mean  rate  of  strain  direction.  However,  in  the  near¬ 
wall  region,  Bradshaw's  turbulent  shear  stress  equation  reduces  to' the 
form,  again  in  the  nomenclature  of  the  present  analysis, 

#  *  (£)2]  <29> 

which  is  identical  with  the  results  of  Eqs.  (25)  and  (26)  derived  above 
for  the  turbulent  shear  contributions.  Arguments  for  the  preference  of 
either  variance  or  invariance  of  the  turbulent  shear  stress  vector  rela¬ 
tive  to  the  mean  rate  of  strain  direction  are  currently  based  on  slender 
evidence.  The  analytical  studies  by  Nash  (Ref.  19)  and  Cooper  (Ref.  20) 
concerning  three-dimensional  incompressible  turbulent  boundary  layers 
indicate  that  the  above-defined  invariant  model  of  turbulence  results  in 
very  acceptable  calculated  boundary -layer  parameters  relative  to  experi¬ 
ment.  Much  more  work,  especially  carefully  controlled  experiments 
involving  three-dimensional  compressible  turbulent  boundary  layers, 
remains  to  be  done  before  the  question  of  variance  or  invariance  can  be 
completely  resolved.  Until  such  time,  the  above  invariant  turbulence 
model  appears  plausible  for  studying  three-dimensional  turbulent 
boundary-layer  flow. 


2.4  MIXING-LENGTH  MODEL 

The  turbulent  shear  stress  in  a  three-dimensional  turbulent  bound¬ 
ary  layer  as  governed  by  Eqs.  (25)  and  (26)  is  treated  herein  by  the  use 
of  a  two-layer  inner-outer  model  using  Prandtl's  mixing-length  hypothesis 
and  a  modification  of  van  Driest 's  analysis  for  the  near-wall  region.  This 
results  in. a  continuous  distribution  of  the  shear  stress  from  the  laminar 
value  at  the  wall,  through  the  fully  turbulent  region,  reaching  zero  at  the 
outer  edge  of  the  boundary  layer.  The  energy  transport  in  a  turbulent 
boundary  layer  is  treated  in  this  work  through  the  incorporation  of  the 
eddy  conductivity,  k,  into  a  turbulent  Prandtl  number,  Pr*.  Exactly 
the  same  empirical  functions  are  used  as  in  the  two-dimensional  flow 
reported  by  Adams  (Ref.  8).  This  is  justified  by  noting  that  the  scalar 
properties  of  a  turbulence  field  are  unlikely  to  be  affected  by  moderate 
three-dimensionality  because  turbulence  is  inherently  three-dimensional, 
in  nature  for  even  so-called  two-dimensional  flows. 

After  Escudier  (Ref.  28),  Patankar  and  Spalding  (Ref.  2)  recommend 
the  following  variation  of  the  mixing  length,  i*,  across  the  turbulent  two- 
dimensional  boundary  layer  which  is  adopted  for  the  three-dimensional 
case  per  the  above  discussion: 
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(30) 


where  the  values  for  the  various  numerical  constants  are  taken  to  be 
k#  =  0.  435  and  X  =  0.  09.  The  value  of  y  at  the  point  where  the  velocity 
in  the  boundary  layer  is  equal  to  0.  99  of  the  velocity  at  the  boundary - 
layer  outer  edge  is  used  to  define  the  distance  y^.  The  above  choices 
follow  Patankar  and  Spalding  and  result  in  good  agreement  with  the 
mixing-length  model  proposed  by  Maise  and  McDonald  (Ref.  29)  for 
compressible  two-dimensional  turbulent  boundary  layers. 


By  analogy  with  Stokes'  solution  for  an  infinite  flat  plate  undergoing 
simple  harmonic  motion  parallel  to  itself  in  an  infinite  fluid,  van  Driest 
(Ref.  30)  concluded  that  in  the  vicinity  of  a  wall  the  total  shear  stress  in 
a  turbulent  two-dimensional  fluid  should  be  of  the  form 


7 


p  t*y2 


i  - 


exp 


(31) 


which  results  in  an  exponential  damping  of  the  turbulent  part  of  the  shear 
stress  as  the  wall  is  approached  and  yields  exactly  the  laminar  shear 
stress  form,  t  =  p(9u/8y),  at  the  wall.  Although  Eq.  ( 31)  was  originally 
developed  for  incompressible  flow,  it  can  be  applied  to  compressible 
flow  by  application  of  the  suggestion  by  Patankar  and  Spalding  (Ref.  2) 
that  the  local  value  of  shear  stress  be  used  instead  of  the  wall  value  as 
originally  recommended  by  van  Driest  (Ref.  30).  Hence,  by  analogy  of 
Eq.  (31)  with  Eqs.  (25)  and  (26),  the  relationship  for  the  three- 
dimensional  near- wall  shear  stress  as  used  in  the  present  analysis  is 


T 


(32) 


where  the  constant  A*  is  taken  to  be  26.  0  following  the  original  van  Driest 
proposal  (Ref.  30).  Note  that  the  damping  term  in  Eq.  (32)  reflects  the 
application  of  the  local  total  shear  stress  as  opposed  to  the  wall  shear 
stress  of  Eq.  (31)  as  discussed  previously. 

Based  on  Eqs.  (26),  (27),  (30),  and  (32),  the  eddy -viscosity  ex¬ 
pression  for  the  inner  region  is 


(33) 
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and  for  the  outer  region 


(34) 


with  the  constants  k#,  A *,  X,  and  y^  defined  previously.  The  constraint 
used  to  define  the  end  of  the  inner  region  and  the  beginning  of  the  outer 
region  is  the  continuity  of  the  eddy  viscosity.  From  the  wall  outward, 
the  expression  for  the  inner  eddy  viscosity  applies  until  =  e0  from 
which  point  the  outer  eddy  viscosity  is  used.  A  schematic  of  this  varia¬ 
tion  in  terms  of  the  mixing  lengths  is  shown  below. 


The  turbulent  Prandtl  number  defined  as 

c  f 

Pr.  =  (35) 

1  K 

is  physically  a  measure  of  the  ratio  of  eddy  viscosity  to  eddy  conduc¬ 
tivity,  that  is,  the  ratio  of  the  turbulent  transport  of  momentum  to  the 
turbulent  transport  of  heat.  Since  the  flow  in  the  outer  region  of  a  tur¬ 
bulent  boundary  layer  shows  some  similarity  to  a  turbulent  wake  flow. 
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one  may  advance  arguments  that  a  realistic  formulation  of  turbulent 
Prandtl  number  requires  a  separate  expression  for  the  inner  and  outer, 
region  just  as  in  the  eddy -viscosity  formalism.  Recent  reviews  by 
Rotta  (Refs.  31  and  32)  and  Meier  and  Rotta  (Ref.  33)  indicate  that  the 
turbulent  Prandtl  number  varies  in  magnitude  from  approximately  0.  8 
near  the  boundary -layer  outer  edge  to  approximately  1.  5  near  the  wall. 
Much  experimental  work  remains  to  be  done  in  defining  the  turbulent 
Prandtl  number  distribution  as  the  reviews  by  Rotta  point  out.  In  the 
present  work  the  turbulent  Prandtl  number  defined  by  Eq.  (35)  is  taken 
to  remain  constant  at  the  value  0.  90  across  the  entire  boundary  layer 
as  recommended  by  Patankar  and  Spalding  (Ref.  2)  for  two-dimensional 
turbulent  boundary  layers.  This  choice  of  turbulent  Prandtl  number  is 
consistent  with  the  analyses  of  Patankar  and  Spalding  (Ref.  2),  Martellucci, 
Rie,  and  Sontowski  (Ref.  6),  Mayne  and  Dyer  (Ref.  7),  and  Adams  (Ref.  8); 
the  analysis  by  Smith  and  Cebeci  (Ref.  1)  assumed  a  value  of  unity. 
Bushnell  and  Beckwith  (Ref.  5),  as  well  as  Hunt,  Bushnell,  and  Beckwith 
(Ref.  17),  employ  a  variable  turbulent  Prandtl  number  model. 


2.5  SUMMARY  OF  GOVERNING  THREE-DIMENSIONAL  TURBULENT 
BOUNDARY-LAYER  EQUATIONS  AND  GAS  MODEL 

Under  the  assumption  of  the  invariant  turbulence  scalar  eddy- 
viscosity  model  of  three-dimensional  turbulence  presented  above,  the 
governing  three-dimensional  turbulent  boundary-layer  equations  for  a 
sharp  cone  at  incidence  flows  under  supersonic  or  hypersonic  conditions, 
as  given  by  Eqs.  (1),  (2),  (3),  (4),  and  (9),  reduce  to  the  following: 


CONTINUITY 


^(plu)  +  ^  (pw)  +  *-(pV r)  =  0 


(36) 


x-MOMENTUM 


.  .  <?u  aw  du  -  y  <Ju 
pu  —  +  —  —  +  P'  -r- 


p(w)2  dr 


dx  r  <30  dy  r  dx 


—  _  j/  +  4-\fa  -  2)- 

dx  Sox  dy  |V  /  dy 


(37) 


0-MOMENTUM 


. .  dw  pw  i?w  -,,<?«  p  uw  dr 

p"5 7*T3?  +  fVW*—X 


-1  db  d  [/  \  divl 


(38) 
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y-momentum 


(39) 


ENERGY 


-  -  r „ 
pucpd7  +  ~  c 


py 

r 


dr 

d<f> 


P'VC 


dl_ 
P  dy 


+ 


2 


(40) 


with  the  scalar  eddy  viscosity  e  defined  by  Eq.  (27)  as 


in  conjunction  with  the  following  mixing-length  model  from  Eqs.  (33) 
and  (34) 

Inner  Region  •»  P*  =■  Dk,y 


Outer  Region  -»  I','  =  Ayp 


Switchover  When 

employing  the  following  empirical  constants 

k*  =  0.435 
A  =  0.090 

D  =  van  Driest  wall  damping  term  evaluated  locally  using  A*  =  26.0 


\i  A, 

yp  =  y-value  where 


[  ( u ) 2  +  ( w)‘] 


‘4 


=  0.99 


[U. 


w  2] 


21*4 


Pq  =  0.90 


and  from  Eq.  (6), 


(41) 
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The  gas  model  used  is  thermally  and  calorically  perfect  air  obeying  the 
equation  of  state  given  by  Eq.  ( 12): 


p  =  pRT 

where 

R  =  1716  ft2/sec2-°R 
Cp  =  6006  ft2/sec2-°R 

y  =  1.40 
Pr  =  0.71 


The  laminar  viscosity  is  given  by  Sutherland's  law,  Eq.  (14),  with  T  in 
units  of  °R: 

_jt _ Tre{+  198.6  /  t  \  * 

^  ret  T+198.S  V'r.f/ 


From  the  sharp  cone  geometry  as  shown  in  Fig.  1, 

r  =  x  sin 

SO  that 


(42) 

(43) 


which  provides  the  necessary  geometry  information  for  use  in  the 
governing  equations  of  motion  (36-40). 


2.6  PROCEDURE  FOR  NUMERICAL  SOLUTION  OF  THE  THREE-DIMENSIONAL 

BOUNDARY-LAYER  EQUATIONS 

For  application  in  the  present  sharp  cone  at  incidence  investiga¬ 
tion,  the  three-dimensional  conical  flow  laminar  boundary -layer  analysis 
as  presented  in  Appendix  B  of  the  report  by  McGowan  and  Davis  (Ref.  34) 
has  been  modified  to  include  the  effects  of  three-dimensional  turbulence 
through  the  use  of  the  scalar  eddy-viscosity  model  discussed  previously. 
The  basic  McGowan  and  Davis  laminar  boundary -layer  treatment  is 
very  similar  to  that  of  Dwyer  (Ref.  35)  and  Boericke  (Ref.  36)  in  that 
the  limiting  forms  of  the  full  three-dimensional  compressible  laminar 
boundary -layer  equations  for  conical  flow  as  originally  derived  by 
Moore  (Ref.  25)  are  solved  using  a  marching  implicit  finite -difference 
technique  for  numerical  integration  of  the  nonlinear  parabolic  partial 
differential  equations  written  in  similarity  variable  form.  The  similar¬ 
ity  variable  transformation  reduces  the  number  of  independent  variables 
from  three  to  two  in  the  transformed  governing  equations  so  that  the 
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problem  becomes  two-dimensional  in  form.  Since  there  are  only  two 
independent  variables  in  this  coordinate  system,  the  implicit  finite- 
difference  techniques  developed  by  Blottner  (Refs.  37  and  38)  can  be 
used  almost  directly  to  solve  the  governing  equations.  The  complete 
formalism  of  this  numerical  approach  is  discussed  in  Chapter  III  of  the 
report  by  McGowan  and  Davis,  to  which  the  reader  is  referred  for 
further  information. 


The  necessary  outer-edge  conditions  for  input  to  the  above- 
described  boundary -layer  analysis  are  determined  based  on  results 
from  an  inviscid  analysis  of  a  sharp  cone  at  incidence,  which  is  dis¬ 
cussed  in  the  following  subsection.  The  procedure  for  specifying  the  in¬ 
viscid  data  necessary  for  input  to  the  McGowan  and  Davis  boundary- 
layer  analysis  is  quite  simple  in  that  only  the  pressure  distribution 
around  the  cone,  along  with  the  velocity  and  density  on  the  windward 
streamline,  must  be  specified.  All  other  inviscid  quantities  are  then 
internally  calculated  using  the  inviscid  compressible  Bernoulli  and 
crossflow  momentum  equations  applied  at  the  cone  surface,  along  with 
the  restriction  that  the  entropy  remain  constant  on  the  surface;  i.  e.  , 
the  cone  surface  is  an  isentropic  surface.  Complete  details  of  this  pro¬ 
cedure  are  given  in  Section  B  of  Chapter  IV  in  the  report  by  McGowan 
and  Davis  (Ref.  34). 


Following  Appendix  B  of  McGowan  and  Davis  (Ref  34),  the  govern¬ 
ing  three-dimensional  turbulent  boundary -layer  equations  (36-40)  are 
transformed  using  similarity  variables  5,  rj,  and  £  similar  to  those  used 
by  Dwyer  (Ref.  35)  and  Boericke  (Ref.  36)  for  three-dimensional  lami¬ 
nar  boundary  layers.  The  definitions  of  C,  rj,  and  £  are  as  follows: 

f  =  /  r2dx  -  f  (x  sin  Sv)2dx  =  -^x3  sin2  8v  (44) 


C  =  <t> 


(45) 

(46) 


where  r  =  x  sin  6V  for  the  sharp  cone  geometry  of  present  interest,  as 
shown  in  Fig.  1.  Introducing  the  above  similarity  variables  into  the 
governing  equations  (36-40)  and  performing  the  standard  transformation 
of  variables  manipulations  yields  the  set  of  equations  (B.  13-B.  16)  in 
Appendix  B  of  McGowan  and  Davis  with  the  following  two  modifications: 


1.  The  laminar  viscosity,  M,  must  be  replaced  by  the 
sum  of  the  laminar  and  turbulent  (eddy)  viscosity 
(M  4  e)  in  the  transformed  §-  and  ?-  momentum  equa¬ 
tions,  as  well  as  in  the  transformed  energy  equation. 
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Furthermore,  in  the  energy  equation  the  laminar  heat 
conductivity  term  (ju/Pr)  must  be  replaced  by  the  sum 
of  the  laminar  and  turbulent  (eddy)  heat  conductivity 
O/Pr)  +  (e/Prt)]. 

2.  The  three-dimensional  turbulent  boundary -layer  flow 
must  be  locally  similar  in  the  sense  of  a  mathematical 
analysis  under  the  constraint  0/8?)  =  0  with  the  eddy 
viscosity  e  evaluated  at  the  local  ?  condition.  The 
applicability  of  this  technique  relies  essentially  on  the 
condition  that  the  external  and  body  flow  properties  vary 
sufficiently  slowly  with  the  x-dependent  variable  ?  de¬ 
fined  by  Eq.  (42).  Experimental  justification  for  the  use 
of  this  assumption  in  the  case  of  three-dimensional  tur¬ 
bulent  boundary-layer  flow  over  a  sharp  cone  at  inci¬ 
dence  in  a  supersonic  stream  is  presented  in  Section  IV 
of  the  present  report. 

Under  the  above  local  similarity  restriction  the  transformed  governing 
boundary-layer  equations  become  mathematically  parabolic  in  the  b, 

£  coordinates  with  ?  as  a  parameter.  The  "history"  of  the  flow  is  con¬ 
tained  only  in  the  e  and  b  dependence  on  ?  and,  hence,  local  similarity 
represents  a  "patching  together"  of  local  solutions. 

The  method  for  numerical  solution  of  the  governing  three- 
dimensional  boundary -layer  equations  in-similarity  £,  b,  £  variables 
follows  the  iterative  implicit  finite-difference  integration  technique  (in¬ 
tegration  in  b-direction  marching  in  £  -direction  windward  to  leeward 
ray)  presented  in  Chapter  III  of  the  report  by  McGowan  and  Davis 
(Ref.  34).  A  variable  b  grid  mesh  is  used  following  Adams  (Ref.  8)  to 
concentrate  grid  points  in  the  near-wall  region  where  the  dependent 
variables  change  most  rapidly  in  a  turbulent  flow.  Digital  computer 
run  times  are  acceptable  for  practical  usage  (approximately  20  minutes, 
including  printout,  to  integrate  180.  0  degrees  around  a  sharp  cone  at 
incidence  using  a  2.50-deg  step  size  ? -integration  increment  with  120  b 
grid  points  across  the  boundary  layer  on  a  CDC  1604-B  digital  computer). 

Experience  with  the  McGowan  and  Davis  digital  computer  code  re¬ 
ported  in  Ref.  34  has  revealed  few  defects,  and  the  present  author 
highly  recommends  its  use.  It  should  be  noted  that  the  main  emphasis 
of  Ref.  34  is  placed  upon  development  and  documentation  of  a  very 
general  three-dimensional  laminar  boundary-layer  analysis  for  general 
body  geometry,  providing  the  inviscid  flow  field  for  the  body  in  question 
is  available  from  some  source. 
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2.7  THREE-DIMENSIONAL  INVISCID  CONICAL  FLOW 

A  recent  investigation  by  Jones  (Ref.  39)  resulted  in  an  accurate 
and  efficient  numerical  integration  procedure  for  solution  of  the  govern¬ 
ing  partial  differential  equations  describing  the  supersonic'or  hyper¬ 
sonic  inviscid  flow  field  around  a  sharp  cone  at  incidence.  Basically, 
Jones'  method  uses  the  condition  of  conicity  to  reduce  the  problem  to  a 
set  of  elliptic  nonlinear  partial  differential  equations  in  two  independent 
variables.  A  transformation  of  coordinates  is  used  to  fix  the  bound¬ 
aries,  one  of  which  is  the  unknown  shock  wave,  between  which  the 
elliptic  equations  are  to  be  satisfied.  This  transformation  also  has  the 
effect  of  including  the  body  shape  in  the  coefficients  of  the  partial  differ¬ 
ential  equations  and  in  the  boundary  conditions,  so  that  the  same  method 
can  be  used  for  general  conical  body  shapes  simply  by  changing  a  few 
program  statements  to  redefine  the  equation  of  the  body.  ■  In  fact,  the 
method  is,  in  many  cases,  only  limited  by  the  crossflow  velocity  ex¬ 
panding  from  subsonic  to  supersonic  conditions  which  changes  the 
mathematical  character  of  the  governing  equations  from  elliptic  to 
hyperbolic,  by  the  entropy  singularity  moving  too  far  away  from  the 
surface,  or  by  the  shock  approaching  very  close  to  the  Mach  wave.  In 
practice  these  restrictions  limit  the  allowable  angle -of -attack  range  to 
al6v  £  1.4  (see  Fig.  1  for  clarification  of  nomenclature). 

The  method  is  efficient  in  computer  time  compared  with  other  fully 
numerical  techniques.  One  solution  takes  from  about  one-half  minute 
to  three  minutes  on  an  IBM  360/50  computer  for  the  circular  cone  at 
incidence,  the  time  increasing  as  the  incidence  increases.  This  is  to 
be  compared  with  a  time  requirement  of  approximately  one-half  hour 
on  an  IBM  360/50  computer  for  the  technique  developed  by  Moretti 
(Ref.  40)  in  which  the  flow  field  solution  is  obtained  by  marching  step 
by  step  downstream  (approximately  400  downstream  steps  are  required) 
until  a  conicity  condition  is  sufficiently  well  satisfied.  Comparison  of 
results  between  the  Jones  and  Moretti  approaches  show  excellent  agree¬ 
ment,  with  the  Jones  digital  computer  code  being  a  factor  of  approxi¬ 
mately  ten  faster  than  the  Moretti  approach  in  solution  time.  An 
analysis  very  similar  to  that  of  Jones  has  recently  been  reported  by 
South  and  Klunker  (Ref.  41),  whereas  Holt  and  Ndefo  (Ref.  42)  have 
developed  a  method  of  integral  relations  approach  to  the  problem.  The 
important  point  to  note  is  that  all  of  the  above -referenced  analyses  re¬ 
port  excellent  agreement  with  experiment  for  sharp  circular  cones  at 
incidence  under  supersonic  and  hypersonic  flow  conditions,  so  that  the 
choice  of  which  analysis  is  indeed  the  best  remains  an  open  question. 

The  present  author's  experience  with  use  of  the  Jones  digital  computer 
code  (Ref.  43)  has  been  most  favorable  from  a  user’s  standpoint. 
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It  should  be  pointed  out  in  conclusion  that  Jones  (Ref.  44)  has 
recently  published  a  very  complete  and  thorough  set  of  tables  for  in- 
viscid  supersonic  and  hypersonic  flow  about  circular  cones  at  incidence 
in  a  perfect  gas,  7  =  1.40,  stream.  These  tables  can  be  used  to  pro¬ 
vide  all  of  the  needed  inviscid  information  for  input  to  the  present 
boundary-layer  analysis. 


SECTION  III 

PRIOR  EXPERIMENTAL  INVESTIGATIONS  OF  THE  SHARP  CONE  AT  INCIDENCE 
UNDER  SUPERSONIC  AND  HYPERSONIC  FLOW  CONDITIONS 


The  primary  set  of  experimental  data  used  for  comparison  with  the 
present  theory  is  taken  from  the  work  of  Rainbird  (Ref.  45)  concerning 
turbulent  boundary -layer  growth  and  separation  on  a  yawed  12.  5 -deg 
semivertex  angle  sharp  cone.  The  investigation  was  conducted  in  the 
Canadian  National  Aeronautical  Establishment  5 -ft  intermittent  blow¬ 
down  (air)  supersonic  wind  tunnel  at  moderate  relative  incidence 
(a! 6V  =  1  to  2)  under  high  Reynolds  number  conditions  with  essen¬ 
tially  zero  heat  transfer  either  to  or  from  the  cone.  Reference  45  pre¬ 
sents  experimentally  determined  surface  pressure  distributions,  sur¬ 
face  flow  angles,  and  detailed  turbulent  boundary -layer  profile  traverses 
at  various  circumferential  locations  around  the  cone.  Full  details  of 
the  experimental  techniques  are  given  in  Ref.  45,  as  well  as  in  an  ear¬ 
lier  paper  by  Rainbird  (Ref.  46)  on  essentially  the  same  subject. 

For  the  present  investigation,  attention  is  restricted  solely  to  the 
free-stream  Mach  Number  1.  80  and  angle-of-incidence  15.  78-deg  con¬ 
dition  of  Rainbird  (Ref.  45).  The  nominal  free-stream  conditions  for 
this  case  are  as  follows  (see  Fig.  1  for  the  sharp  cone  geometry  and 
nomenclature): 

IV1  =  1.80 

oo 

p*.  =  626.40  lbt/ft2 

T  =  321.60°R 

OO 

ReK  L  =  2.56  x  107 

L  =  41.58  in. 

The  cone  surface  temperature  is  taken  as  equal  to  the  free-stream  stag¬ 
nation  temperature  (530°R).  Because  of  the  conical  nature  of  the  flow 
field  for  a  sharp  cone  at  incidence  (discussed  in  detail  in  Section  4.  1  of 
this  report),  all  boundary -layer  surveys  at  various  circumferential 
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locations  around  the  cone  were  conducted  at  one  axial  location  along  the 
cone,  namely,  x/L  =  0.85.  As  stated  by  Rainbird  (Ref.  45),  boundary- 
layer  transition  takes  place  quite  close  to  the  sharp  cone  apex  (x/L  <  0.1) 
because  of  the  high  stream  turbulence  level  resulting  from  noise  gen¬ 
erated  by  the  blowdown  wind  tunnel  control  valve. 

The  other  set  of  experimental  data  used  for  comparison  with  the 
present  theory  is  taken  from  the  study  by  Martellucci  and  Neff  (Ref.  47) 
concerning  the  effects  of  asymmetric  boundary -layer  transition  on  re¬ 
entry  vehicle  characteristics  under  hypersonic  flow  conditions.  The 
experimental  phase  of  this  investigation  was  conducted  in  the  Hypersonic 
Wind  Tunnel  (B)  of  the  von  Karman  Gas  Dynamics  Facility  (VKF),  Arnold 
Engineering  Development  Center  (AEDC),  on  a  7.  20-deg  semivertex 
angle  sharp  cone  at  1.  0-,  2.0-,  and  4.  0-deg  angle  of  attack.  The  model, 
which  was  fabricated  of  stainless  steel,  had  a  nominal  wall  thickness  of 
0.  050  in.  and  was  instrumented  with  approximately  100  thermocouples 
located  on  five  conical  rays  (0=  0,  45,  90,  135,  and  180  deg)  to  record 
the  aerodynamic  heat-transfer  rate  via  the  thin-skin  technique.  Tunnel 
free-stream  conditions  are  as  follows  (see  Fig.  1  for  the  sharp  cone 
geometry  and  nomenclature): 

=  8.0 

=  12.69  Ibf/ft2 

=  96.84°R 

Re«.lL  =  1-38  x  107 
L  =  43.89  in. 

which  are  based  upon  nominal  tunnel  stagnation  conditions  of 

p0  -  860  Ibf/in." 

T0  =  1340°R 

The  cone  surface  temperature  is  taken  to  be  a  constant  535.  0°R,  which 
results  in  a  constant  wall  to  stagnation  temperature  ratio,  TW/TQ,  of 
approximately  0.  4.  Figure  2  presents  the  spatial  distribution  of  the 
end  of  boundary -layer  transition  with  respect  to  angle  of  attack  as  re¬ 
ported  in  Fig.  3  of  Ref.  47  for  the  above  tunnel  free-stream  conditions. 
Here  the  end  of  transition  is  defined  to  be  the  point  where  the  aero¬ 
dynamic  heat -transfer  rate  in  the  transition  front  region  is  a  maximum. 
The  important  point  to  note  from  Fig.  2  is  that  the  three-dimensional 
boundary  layer  will  be  in  a  fully  turbulent  state  for  x/L  >  0.  5  for  all 
three  angles  of  attack  under  the  present  high  Reynolds  number  free- 
stream  condition. 
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4.1  JUSTIFICATION  OF  LOCALLY  SIMILAR  TURBULENT 
BOUNDARY-LAYER  ANALYSIS 

The  assumption  made  in  Section  II  of  a  locally  similar  boundary- 
layer  analysis  with  the  eddy  viscosity  e  evaluated  at  the  local  £  condi¬ 
tion  appears  very  questionable  for  application  to  general  three- 
dimensional  turbulent  flows  because  of  the  failure  to  include  details  of 
the  "upstream  history. "  However,  for  the  special  case  of  a  sharp  cone 
at  incidence,  in  a  supersonic  or  hypersonic  stream  where  the  boundary, 
layer  is  in  a  state  of  fully  developed  turbulent  flow,  i.  e. ,  far  down¬ 
stream  of  transition  with  a  constant  wall  temperature,  experimental 
measurements  reported  by  Rainbird  (Refs.  45  and  46)  establish  that  the 
flow  field,  even  with  separation  present,  is  essentially  conical  and  sym¬ 
metrical,  thus  permitting  all  detailed  measurements  to  be  made  at  one 
lengthwise  station.  The  evidence  in  support  of  this  finding  is  as  follows 
(taken  from  Ref.  45): 

a.  Overall  force  and  moment  measurements  show  zero 
side  force  and  yawing  moment  and  give  a  fixed 
center -of-pressure  position  at  0.  682L  that  is  in  ex¬ 
cellent  agreement  with  the  theoretical  conical  flow 
value  of  (2/3)L/cos  6V. 

b.  Integration  of  circumferential  pressure  distributions 
to  give  local  normal-force  coefficients  shows  good 
agreement  with  overall  balance  measurements. 

c.  Measurements  of  surface  pressure  distributions  along 
generators  of  the  cone  show  pressures  constant  ex¬ 
cept  for  some  extreme  angle -of -attack  conditions 
where  a  forward-propagating  base  effect  is' present. 

d.  Flow  visualization  using  the  oil-dot  technique  gives 
values  of  surface  flow  angle  us,  i.  e.  ,  the  direction 
of  surface  shear  stress,  as  well  as  primary  separa¬ 
tion  position  which  are  independent  of  distance  from 
the  cone  apex,  x/L,  within  a  measuring  accuracy  of 
about  1.  5  deg. 

Because  of  the  importance  of  the  invariance  of  the  surface  flow 
angle  with  lengthwise  location  at  a  given  circumferential  location  in  the 
present  theoretical  analysis,  a  comparison  is  given  in  Fig.  3  (taken 
from  Ref.  45)  of  surface  flow  angle  measurements  at  various  x/L  sta¬ 
tions  up  to  separation  for  the  rather  severe  condition  of  a! 6W  =  2.  See 
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Fig.  1  for  the  definition  of  the  surface  flow  angle  relative  to  the  conical 
geometry  of  present  interest.  The  results  show  that  surface  flow  angle 
is  essentially  independent  of  distance  from  the  cone  apex,  which  means 
that  under  such  a  flow  condition  {conical  inviscid  and  fully  turbulent 
boundary  layer)  a  locally  similar  turbulent  boundary-layer  analysis 
which  neglects  "upstream  history"  may  be  a  plausible  assumption.  One 
purpose  of  the  present  report  is  to  assess  if  indeed  such  is  a  reasonable 
approach  based  on  comparison  of  theory  relative  to  experimental  mea¬ 
surements  taken  under  flow  conditions  which  satisfy  the  restrictions 
imposed  above  (conical  inviscid  flow  field  with  a  fully  turbulent  bound¬ 
ary  layer  under  constant  wall  temperature  conditions). 


4.2  PRESENTATION  OF  PRESENT  RESULTS 
4.2.1  Supersonic  Flow 

Turning  now  to  representative  results  from  the  present  investiga¬ 
tion,  one  sees  in  Fig.  4  a  comparison  of  the  calculated  surface  pressure 
distribution  around  the  sharp  cone  based  on  the  Jones  analysis  (Refs.  39 
and  43)  relative  to  the  experimental  measurements  of  Rainbird  (Ref.  45). 
As  is  clearly  shown  in  Fig.  4,  the  agreement  is  excellent  over  the  en¬ 
tire  cone.  Figure  5  presents  the  corresponding  calculated  inviscid  flow 
parameters  (streamwise  and  crossflow  velocities,  as  well  as  static 
temperature)  on  the  cone  surface.  It  should  be  noted  that  these  surface 
values  are  the  so-called  isentropic  surface  values  (see  Ref.  39  for 
clarification  of  this  terminology). 

With  respect  to  the  above,  it  should  be  noted  that  Rainbird  (Ref.  45) 
experimentally  observed  turbulent  boundary -layer  separation  to  occur 
at  approximately  159  deg  around  the  cone  for  the  present  flow  condition 
and  angle  of  incidence  (MD  =  1.  80,  Re^  ^  =  2.  56  x  10?,  a  =  15.  78  deg). 
As  discussed  by  Rainbird  in  Ref.  45,  ttie  development  of  flow  separa¬ 
tion  about  sharp  cones  as  the  incidence  angle  is  increased  is  a  gradual, 
progressive,  steady,  and  essentially  conical  process  involving  the 
formation  of  symmetrical  lobes  of  vortical  fluid  which  develop,  into 
vortices  arjd  which  remain  comparatively  close  to  the  cone  surface  on 
either  side  of,  and  near,  the  leeward  generator.  For  the  angle  of  inci¬ 
dence  of  present  interest  ( a  =  15.  78  deg)  Rainbird  observed  the  forma¬ 
tion  of  two  symmetrically  disposed  lobes  of  vortical  fluid  on  either  side 
of  the  leeward  generator.  At  a  higher  angle  of  incidence  (a  =  22.  75  deg) 
these  lobes  of  vortical  fluid  roll  up  to  form  a  pair  of  symmetrically  dis¬ 
posed  vortices  close  to  the  cone  surface  which,  in  turn,  result  in  the 
formation  of  internal  shock  waves  with  their  attendant  local  increases 
in  pressure.  Since  there  are  no  vortices  present  in  the  separated  flow 
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field  of  current  interest,  the  influence  of  separation  on  the  external  in- 
viscid  flow  is  small,  which  is  reflected  in  the  excellent  agreement 
shown  in  Fig.  4  between  inviscid  theory  and  experiment. 

Figures  6  and  7  present  the  calculated  mean  velocity  and  static 
temperature  profiles  across  the  turbulent  boundary  layer  at  the  location 
x/L  =  0.  85  on  the  most  windward  ray  (tf>  =  0  deg)  of  the  sharp  cone.  The 
calculated  profiles  are  generally  in  good  agreement  with  the  measured 
profiles  by  Rainbird,  which  reveals  the  validity  of  the  presently  pro¬ 
posed  three-dimensional  eddy-viscosity  model  for  windward  ray  applica¬ 
tions.  Shown  in  Fig.  7  are  comparisons  of  the  present  windward  ray 
profiles  relative  to  calculated  results  from  Fig.  16b  of  the  recent  wind¬ 
ward  plane  of  symmetry  turbulent  boundary-layer  analysis  by  Adams 
(Ref.  48).  The  basic  approach  of  Ref.  48  involved  formulation  and 
application  of  a  laminar,  transitional,  and  turbulent  boundary -layer 
analysis  for  the  windward  streamline  of  a  sharp  cone  at  incidence  in  a 
supersonic  or  hypersonic  flow.  The  governing  nonsimilar  boundary- 
layer  equations  in  the  windward  plane  of  symmetry  were  numerically 
integrated  on  a  digital  computer  using  an  implicit  finite -difference  tech¬ 
nique  which  marched  along  the  windward  ray  starting  at  the  apex  of  the 
cone  with  a  laminar  similar  solution.  The  same  two-layer  (inner-outer) 
eddy  viscosity-mixing  length  model  of  turbulence  was  used  for  calcula¬ 
tion  of  the  windward  ray  turbulent  boundary  layer  as  in  the  present  work. 
The  transition  zone  was  treated  through  an  eddy  viscosity-intermittency 
factor  approach.  Inviscid  edge  conditions  along  the  windward  ray  were 
obtained  from  the  same  Jones  digital  computer  code  (Refs.  39  and  43) 
used  in  the  present  work.  The  excellent  agreement  shown  in  Figs.  7a 
and  b  between  the  nonsimilar  analysis  of  Ref.  48  and  the  present 
locally  similar  analysis  offers  further  analytical  justification  for  the  ■ 
applicability  of  the  locally  similar  type  analysis  for  sharp  cone  at  inci¬ 
dence  flows  with  a  turbulent  boundary  layer.  It  should  be  mentioned  in 
this  connection  that  Schmidt,.  Boldman,  and  Todd  (Ref.  49)  have  re¬ 
cently  reported  a  locally  similar  turbulent  boundary-layer  analysis  for 
axisymmetric  nozzle  flow  applications.  They  point  out  that  through 
this  approach  the  nozzle  exit  boundary  layer  can  be'  calculated  directly 
without  the  requirement  of  the  usual  step-by-step  marching  calculation 
method  of  the  nonsimilar  boundary  layer.  The  same  statement  applies 
to  the  present  sharp  cone  at  incidence  boundary  layer  where,  as  in  the 
current  work,  the  analysis  can  be  performed  for  the  body  station  loca¬ 
tion  of  interest  without  having  to  perform  the  step-by-step  marching 
calculation  of  the  complete  three-dimensional  nonsimilar  boundary 
layer. 

Also  shown  in  Figs.  7a  and  b  are  rough  indications  of  the  physical 
boundaries  for  the  inner  and  outer  regions  of  the  turbulent  boundary 
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layer,  as  discussed  by  Bradshaw  (Ref.  50)  in  his  recent  review.  The 
location  of  the  inner-outer  boundary  at  y/6t«  0.  2  can  easily  be  derived 
from  Eq.  (30)  by  equating  the  inner  and  outer  mixing  lengths  to  yield 


_y 

n 


_A_ 

k* 


0.090 

0.435 


0.206 


(47) 


where  y^  is  the  value  of  y  at  the  point  where  the  velocity  in  the  boundary 
layer  is  equal  to  0.  99  of  the  velocity  at  the  boundary-layer  outer  edge, 
and,  hence,  y#  «  6^  where  6*  is  the  boundary-layer  thickness  defined  as 
the  normal  distance  from  the  surface  where  the  boundary-layer  velocity 
equals  0.  995  of  the  outer-edge  velocity.  Further,  observe  from 
Figs.  7a  and  b  that  the  experimentally  determined  turbulent  boundary- 
layer  thickness  is  slightly  larger  than  predicted  by  the  present  theory. 


Using  the  implicit  finite -difference  integration  technique  to  obtain 
the  solution  around  the  cone  at  the  body  location  x/L  =  0.  85  yields  the 
calculated  profiles  shown  in  Figs.  8a,  b,  and  c  for  the  angular  locations 
$  =  45.0,  90.0,  and  135.0  deg,  respectively.  As  can  be  seen  from 
these  figures,  agreement  between  the  calculated  profiles  and  the  experi¬ 
mental  profiles  are  in  good  agreement  for  the  $  =45.  0-deg  case,  differ 
somewhat  in  the  near- wall  region  for  the  0  =  90.  0-deg  case,  and  differ 
somewhat  in  character  across  the  entire  profile  at  <j)  =  135.  0  deg. 

Figure  9  presents  the  streamline  direction  within  the  boundary  layer 
which  shows  good  agreement  between  calculated  and  measured  values 
for  the  <f>  =  45.  0-deg  case  and  progressive  disagreement  as  the  (f>  -angle 
is  increased.  This  behavior  can  be  partially  traced  to  the  use  of  the 
isentropic  surface  values  of  the  inviscid  flow  quantities  as  the  boundary- 
layer  outer-edge  conditions.  As  discussed  by  Rainbird  (Ref.  45)  in  the 
concluding  paragraph  of  his  paper,  it  is  perhaps  more  appropriate  to 
use  "near"  surface  conditions  (rather  than  isentropic  surface  conditions) 
as  the  external  flow  for  bound  ary -layer  calculations.  More  work  re¬ 
mains  to  be  done  in  this  connection. 


Distributions  of  the  surface  flow  angle,  ug,  and  the  external  flow 
angle,  ue,  relative  to  experimental  measurements  are  presented  in 
Fig.  10.  The  condition  that  us  =  0  on  a  conical  surface  is  used  as  a 
criterion  for  boundary -layer  separation  following  Rainbird  (Refs.  45 
and  46).  As  can  be  seen  from  Fig.  10,  the  present  three-dimensional 
turbulent  boundary-layer  analysis  predicts  separation  to  occur  some¬ 
where  between  $  =  162.  5  and  165.0  deg,  whereas  Rainbird  (Ref.  45) 
experimentally  observed  separation  at  <t>  »  159  deg.  Further  note  that 
the  magnitude  of  the  crossflow  influence  on  the  turbulent  bound  ary  -layer 
turning  is  very  small,  e.  g. ,  us  -  ue  »  7  deg  at  <j>  -  90.  0  deg.  .  Also 
shown  in  Fig.  10  is  the  calculated  surface  flow  angle  distribution  for  a 
laminar  boundary  layer  under  the  same  flow  conditions.  Much  larger 
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crossflow  influence  on  the  laminar  bound  ary -layer  turning  is  observed, 
i.  e. ,  us  -  ue  »  26  deg  at  <j>  =  90.  0  deg.  Laminar  boundary-layer  sepa¬ 
ration  is  predicted  to  occur  much  earlier  than  for  the  turbulent  case, 
somewhere  between  $  =  130.  0  and  132.  5  deg.  No  attempt  has  been 
made  in  the  present  study  to  attempt  more  accurate  location  of  the 
calculated  separation  location  by  use  of  a  very  small  <j>  integration  in¬ 
crement  near  separation;  all  of  the  present  calculations  employed  a 
constant  $  integration  increment  of  2.  50  deg.  The  above  results  con¬ 
firm  quite  clearly  the  statement  by  Vaglio-Laurin  (Ref.  11)  that  "due 
to  the  larger  shearing  stress,  smaller  three-dimensional  effects  can 
be  expected  for  turbulent  layers  as  compared  with  laminar  layers  sub¬ 
ject  to  the  same  boundary  conditions.  " 

The  calculated  eddy- viscosity,  e,  distributions  across  the  three- 
dimensional  turbulent  boundary  layer  at  various  angular  locations 
around  the  cone  are  presented  in  Fig.  11.  As  is  apparent  the  eddy 
viscosity  reaches  its  maximum  value  (e  max  «  200  to  300  p)  in  the 
outer  region  of  the  boundary  layer  with  e  »  p  even  in  regions  near  the 
wall;  i.  e.  ,  e  «  lOp  at  y  »  0.  0025  in.  It  should  be  noted  that  the  laminar 
viscosity,  P,  in  the  above  is  evaluated  at  the  same  local  conditions  as 
the  corresponding  eddy  viscosity. 

To  illustrate  the  differences  in  three-dimensional  laminar  and 
turbulent  boundary-layer  structure  at  a  common  body  location.  Figs.  12a, 
b,  and  c  present  calculated  boundary -layer  profiles  at  the  circumferen¬ 
tial  location  $  =  90.  0  deg  for  the  present  sharp  cone  at  incidence  flow. 

The  differences  in  the  laminar  and  the  turbulent  profiles  are  apparent. 

It  is  interesting  to  note  that  the  turbulent  boundary -layer  thickness,  6^-, 
is  approximately  a  factor  of  ten  greater  than  the  laminar  value,  6%,  for 
this  particular  flow  condition. 

4.2.2  Hypersonic  Flow 

Turning  now  to  the  hypersonic  flow  past  a  7.  20 -deg  semivertex 
angle  sharp  cone  at  1.0-,  2.0-,  and  4.0-deg  angles  of  attack  under 
AEDC-VKF  Tunnel  B  (Mach  8)  conditions,  one  sees  in  Fig.  13  the  calcu¬ 
lated  inviscid  flow  parameters  around  the  (isentropic)  cone  surface  based 
on  the  Jones  analysis  (Refs.  39  and  43).  As  is  obvious  from  Fig.  13, 
the  inviscid  parameters  most  sensitive  to  increasing  angle  of  attack  are 
the  surface  pressure,  temperature,  and  crossflow  velocity;  the  stream- 
wise  velocity  component  remains  almost  constant  at  a  value  slightly 
less  than  the  free-stream  velocity  over  the  entire  cone.  The  small 
magnitude  of  the  inviscid  crossflow  velocity  at  even  4 -deg  angle  of 
attack  indicates  that  the  inviscid  crossflow  turning,  i.  e. ,  the  inviscid 
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streamline  direction,  (Je,  at  the  outer  edge  of  the  boundary  layer,  will 
be  small.  More  will  follow  on  this  point  later. 

One  of  the  primary  objectives  of  the  experimental  investigation  re¬ 
ported  by  Martellucci  and  Neff  (Ref.  47)  was  the  detection  of  boundary- 
layer  transition  on  cones  at  angle  of  attack  via  surface  heat -transfer 
measurements,  as  discussed  in  Section  III.  To  establish  the  basic 
validity  of  the  McGowan  and  Davis  (Ref.  34)  analysis  of  the  three- 
dimensional  laminar  boundary  layer  on  a  sharp  cone  at  incidence  under 
hypersonic  conditions.  Fig.  14  presents  a  comparison  of  calculated 
laminar  boundary-layer  heat-transfer  rate  (as  reflected  through  the 
Stanton  number,  St^  based  on  free-stream  conditions)  relative  to 
experimental  measurements  for  a  4.0-deg  angle -of -attack  laminar  flow 
condition  reported  in  Ref.  47.  Note  that  the  heat  transfer  is  presented 
in  laminar  boundary-layer  similarity  format;  i.  e. ,  the  Stanton  number 
is  multiplied  by  the  square  root  of  the  normalized  surface  distance  from 
the  apex.-  Free-stream  conditions  for  this  laminar  case  are  presented  in 
the  figure.  In  general,  the  agreement  between  the  McGowan  and  Davis 
analytical  analysis  and  the  experimental  measurements  of  Ref.  47  is 
excellent  for  this  laminar  flow  condition,  which  illustrates  the  basic 
validity  of  the  McGowan  and  Davis  analysis  for  laminar  boundary  layers 
on  sharp  cones  at  incidence  under  hypersonic  conditions.  It  should  be 
noted  here  that  Boericke  (Ref.  36)  has  also  reported  good  agreement 
with  the  experimental  measurements  of  Ref.  47  for  a  sharp  cone  at 
incidence  under  laminar  boundary-layer  hypersonic  conditions. 

The  end  of  transition  locations  presented  previously  in  Fig.  2 
show  that  the  boundary  layer  will  be  in  a  fully  turbulent  state  for 
x/L  >  0.  5  at  all  three  angle -of -attack  conditions  under  the  high 
Reynolds  number  free-stream  conditions  discussed  in  Section  III.  The 
calculated  turbulent  heat-transfer-rate  distributions  from  the  present 
three-dimensional  turbulent  boundary-layer  analysis  (as  again  reflected 
through  the  Stanton  number  based  on  free-stream  conditions)  for 
x/L,  >  0.  5  and  all  three  angles  of  attack  (a  =  1.  0,  2.  0,  and  4.  0  deg)  at 
various  circumferential  locations  around  the  cone  are  presented  in 
Fig.  15  relative  to  the  experimental  heat-transfer-rate  measurements 
of  Ref.  47.  The  present  calculations  were  performed  at  the  x/L  loca¬ 
tions  indicated  by  the  x-marks  on  the  figure,  and  the  solid  line  was 
faired  between  the  marked  calculated  values.  In  general,  good  agree¬ 
ment  between  the  present  analytical  analysis  and  experiment  is  observed 
at  all  thfree  angle-of-attack  conditions  with  the  circumferential  distribu¬ 
tion  also  well  defined. 
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It  is  well-known  that  in  the  free -stream  Reynolds  number  range 
5  x  10^  <  Re^L  <  10^  the  turbulent  heat-transfer  rate  to  a  flat  plate  in 
incompressible  flow  can  be  expressed  in  the  form  (see,  e.  g. ,  Ref.  51, 
pp.  492-499,  536-539) 


St 


(Re„  J 


-0.2 


(48) 


which  means  physically  that  for  two  x-locations  along  the  plate,  say 
locations  xj  and  X2,  the  corresponding  turbulent  heat-transfer  rates  at 
these  two  locations  are  related  according  to 

<StJ.,  _ 

1S'-’>2  Vi)  (49) 

As  discussed  in  Ref.  51,  Eq.  (49)  may  be  carried  over  reasonably  well 
to  turbulent  compressible  flows,  providing  the  wall  temperature  remains 
constant  and  the  external  flow  at  the  outer  edge  of  the  boundary  layer 
does  not  vary  along  the  body.  Since  these  criteria  are  met  .in  the  pres¬ 
ent  sharp  cone  at  incidence  flow  for  the  case  of  the  streamwise  directed 
velocity  component,  it  is  of  interest  to  note  that  the  scaling  criterion 
suggested  by  Eq.  (49)  works  very  well  when  applied  to  the  results  of 
Fig.  15  in  the  interval  0.  5  <  x/L  <  1.  0  for  all  three  angle -of- attack 
conditions.  This  means  that  for  practical  applications  the  present 
analytical  calculation  of  the  three-dimensional  turbulent  boundary  layer 
need  only  be  performed  at  one  x/L  location,  say  the  x/L  =  0.  80  location, 
as  depicted  in  Fig.  16  for  the  present  flow  condition,  and  then  these  re¬ 
sults  for  the  one  station  scaled  through  the  use  of  Eq.  (49)  to  other  x/L 
locations  at  a  given  circumferential  location.  An  illustration  of  the 
applicability  of  this  approximate  technique  is  given  in  Table  I  for  the 
a  =  4.  0-deg  flow  condition  of  present  interest  with  the  x/L  =  0.8  station 
of  Fig.  16  used  as  the  reference  station  for  the  scaling  application. 
Examination  of  Table  I  reveals  that  for  the  present  body,  angle  of  attack, 
and  flow  condition  the  maximum  error  in  the  use  of  Eq.  (49)  relative  to 
the  present  three-dimensional  turbulent  boundary-layer  analysis 
(denoted  as  "exact"  in  Table  I)  is  only  approximately  2  percent,  which 
is  for  the  x/L  =  0.  50  station  on  the  leeward  ray  ($  =  180.  0  deg). 

Through  the  use  of  this  approximate  scaling  technique,  rapid  calculations 
of  three-dimensional  turbulent  heat-transfer  rates  on  sharp  cones  at 
incidence  can  be  performed  for  use  in  design  tradeoff  studies,  ablation 
heat  shield  design,  or  evaluation  of  experimental  data  with  a  minimum 
requirement  of  digital  computer  time  (one  numerical  calculation  of  the 
three-dimensional  turbulent  boundary  layer  using  the  present  analytical 
analysis).  Needless  to  say  the  applicability  of  the  technique  should  be 
first  established  for  the  flow  conditions  and  geometry  of  interest  in  a 
new  problem,  since  Eq.  (49)  is  not  universally  applicable  for  all  values 
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of  the  free-stream  Reynolds  number  (see  Chapters  XIX  and  XXI  of 
Ref.  51  for  further  discussion  of  this  point). 

Since  the  favorable  agreement  between  results  from  the  present 
theory  and  experiment  indicates  that  the  present  analytical  analysis  is 
indeed  applicable  to  calculation  of  turbulent  heat  transfer  on  a  sharp 
cone  at  incidence  in  hypersonic  flow  under  cold  wall  conditions,  it  is 
now  in  order  to  examine  certain  details  of  the  flow  field.  Figure  17 
presents  the  calculated  windward  ray  ( <j>  ~  0.  0  deg)  turbulent  boundary- 
layer  profiles  at  the  body  station  x/L  =  0.  80  for  the  4.  0-deg  angle-of- 
attack  condition.  Also  included  in  Fig.  17  are  the  calculated  profiles 
based  on  the  nonsimilar,  windward  plane  of  symmetry,  turbulent 
boundary-layer  analysis  by  Adams  (Ref.  48).  The  same  comments 
made  previously  with  respect  to  discussion  of  Figs.  7a  and  b  apply 
equally  as  well  to  Fig.  17.  Further  note  from  Fig.  17  that  the  present 
locally  similar  analysis  and  the  nonsimilar  analysis  of  Ref.  48  are  in 
excellent  agreement  in  the  inner  region  of  the  turbulent  boundary  layer. 
This  is  because  of  the  "local  equilibrium"  and  "independence"  in  the 
inner  layer  of  a  turbulent  boundary  layer,  as  discussed  by  Bradshaw 
(Ref.  50).  A  slight  discrepancy  between  the  two  analyses  is  observed 
in  the  outer  region  of  the  turbulent  boundary  layer,  especially  for  the 
static  temperature.  This  is  because  it  is  the  outer  layer  of  the  turbu¬ 
lent  boundary  layer  where  the  "memory"  of  the  turbulence  becomes 
important  (see  Ref.  50  for  further  discussion  of  this  point). 

The  calculated  turbulent  boundary-layer  parameters  (heat -transfer 
rate,  skin  friction,  and  displacement  thickness)  along  the  windward  ray 
of  the  present  sharp  cone  at  4.  0-deg  angle  of  attack  are  shown  in  Fig.  18. 
The  excellent  agreement  in  heat-transfer  rate  to  the  wall  (through  the 
Stanton  number  St,,,)  and  wall  shear  stress  (through  the  skin-friction 
coefficient  Cf^)  between  the  present  locally  similar  type  of  analysis  and 
the  nonsimilar  analysis  of  Ref.  48  is  directly  due  to  the  excellent  agree¬ 
ment  in  mean  velocity  and  static  temperature  profiles  in  the  inner 
region  of  the  turbulent  boundary  layer,  as  shown  in  the  previous  figure. 
The  small  difference  in  magnitude  of  the  displacement  thickness  ^2-D 

between  the  two  analyses  is  due  to  the  difference  in  static  temperature 
profiles  across  the  outer  region  of  the  turbulent  boundary  layer,  as 
shown  previously  in  Fig.  17.  It  should  be  noted  that  the  displacement 
thickness  definition  used  in  the  present  work  (see  the  NOMENCLATURE) 
is  for  a  two-dimensional  boundary  layer;  see  Adams  (Ref.  48,  p.  22) 
for  a  discussion  of  the  differences  between  two-  and  three-dimensional 
displacement  thickness  definitions  for  the  windward  ray  of  a  sharp  cone 
at  incidence. 
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The  calculated  turbulent  boundary -layer  profiles  for  the  leeward 
ray  (<j>  =  180.  0  deg)  at  the  body  station  x/L  =  0.  80  and  the  4.  0-deg 
angle -of -attack  condition  are  given  in  Fig.  19.  As  originally  shown  by 
Moore  (Ref.  52)  and  recently  discussed  in  great  detail  by  Murdock 
(Ref.  53),  the  laminar  boundary -layer  equations  for  the  leeward  plane 
of  symmetry  do  not  have  a  mathematically  unique  solution  because  of  a 
defect  in  the  basic  boundary-layer  model  on  the  lee  ray  of  a  sharp  cone 
at  incidence.  Furthermore,  Boericke  (Ref.  36)  shows  by  integration  of 
the  three-dimensional  laminar  boundary-layer  equations  around  a  sharp 
cone  at  incidence  from  the  windward  to  the  leeward  ray  that  smooth 
solutions  do  not  exist  in  the  leeward  plane  of  symmetry;  the  same  type 
results  and  conclusions  have  also  been  obtained  by  McGowan  and  Davis 
(Ref.  34,  pp.  32-33).  As  can  be  seen  by  a  careful  study  of  the  analysis 
by  Murdock  (Ref.  53),  the  problem  with  the  leeward  ray  plane  of 
symmetry  bound  ary -layer  equations  can  be  traced  to  the  convective 
terms  in  the  crossflow  momentum  equation  and,  hence,  the  above  dis¬ 
cussion  of  the  laminar  equations  also  applies  to  the  turbulent  case  of 
present  interest,  since  the  only  modifications  for  turbulence  involve 
inclusion  of  the  turbulent  shear  stress  and  conductivity  terms.  Indeed, 
it  has  been  found  in  the  present  work  that  numerical  solutions  could  not 
be  obtained  by  integrating  the  three-dimensional  turbulent  boundary- 
layer  equations  around  a  sharp  cone  at  incidence  from  the  windward  to 
the  leeward  ray— a  converged  solution  in  the  leeward  plane  of  symmetry 
could  not  be  generated.  In  view  of  these  difficulties,  the  profiles  pre¬ 
sented  in  Fig.  19  are  the  results  of  extrapolating  the  numerical  solu¬ 
tions  at  the  two  previous  0  -stations  off  the  leeward  ray  into  the  leeward 
ray  using  a  second-order  extrapolation  process.  While  this  procedure 
is  certainly  not  rigorously  justifiable,  the  favorable  results  presented 
previously  in  Fig.  15  with  respect  to  lee  ray,  turbulent  heat-transfer 
calculations  suggest  that  the  present  extrapolation  technique  yields 
reasonable  turbulent  bound  ary -layer  profiles.  What  is  needed  here  for 
completeness  are  detailed,  lee  ray,  turbulent  boundary -layer  profile 
measurements  for  comparison  with  the  present  analytical  approach. 

Calculated,  three-dimensional,  turbulent  boundary-layer  profiles 
at  three  0  -locations  (0  =  45.0,  90.0,  and  135.0  deg)  are  presented  in 
Fig.  20  for  the  4.  0-deg  angle -of -attack  condition  at  the  body  location 
x/L  =  0.  8.  Note  the  increasing  crossflow  influence  on  the  profiles  as 
the  0  -angle  is  increased,  especially  with  respect  to  the  overshoot  in 
the  crossflow  velocity  profiles  of  Fig.  20b.  Most  of  the  hypersonic 
three-dimensional  turbulent  boundary-layer  turning  relative  to  the 
inviscid  flow,  occurs  in  the  outer  region  (y  >  0.  26t)  of  the  layer,  as  can 
be  seen  from  Fig.  20d.  For  the  0  =  90. 0-deg  case  in  Fig.  20d  with 
6j./L  *8  x  10  "3,  we  »  4  deg,  and  us  *  10  deg  at  y/L  =  10 u  *  8  deg  at 
y/L  =  1.  6  x  IQS,  which  is  the  location  where  y/L  *0.  26t/L.  The 
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three-dimensional  eddy-viscosity  profiles  shown  in  Fig.  20 e  are  very 
similar  in  appearance  to  the  zero-angle-of-attack  eddy-viscosity  pro¬ 
files  presented  in  Fig.  13  of  the  hypersonic  sharp  cone  work  by  Adams 
(Ref.  8).  Further  note  the  somewhat  dissimilar  character  between  the 
eddy -viscosity  profiles  in  Fig.  11  for  the  supersonic  hot  wall  flow  con¬ 
dition  and  the  Fig.  20e  cold  wall  hypersonic  flow  profiles.  This  differ¬ 
ence  in  character  is  directly  related  to  the  differences  in  static  temper¬ 
ature  profiles  between  the  two  flows,  as  can  be  seen  by  comparison  of 
Figs.  8a,  b,  and  c  with  Fig.  20c,  which,  in  turn,  are  related  to  the 
basic  differences  between  a  hot  wall  supersonic  flow  and  a  cold  wall 
hypersonic  flow  over  a  sharp  cone  at  incidence. 

Returning  to  the  subject  of  three-dimensional  boundary -layer  turn¬ 
ing  relative  to  the  inviscid  flow  direction  on  a  sharp  cone  at  incidence 
under  hypersonic  conditions.  Fig.  21  shows  calculated  surface  and 
external  (inviscid)  flow  angles  for  both  laminar  and  turbulent  boundary 
layers  at  the  body  location  x/L  =  0.  8  for  the  4. 0-deg  angle -of- attack 
condition  of  present  interest.  At  the  #  =  110. 0-deg  circumferential 
location  where  maximum  turning  occurs.  Fig.  21  shows  that  us  «  25  deg 
for  the  three-dimensional  laminar  boundary  layer,  whereas  us  «  11  deg 
for  the  three-dimensional  turbulent  boundary  layer;  ue  *  4  deg  for  the 
inviscid  external  flow  at  this  location.  Hence,  the  maximum  laminar 
boundary -layer  turning  angle  is  approximately  a  factor  of  six  greater 
than  the  inviscid  flow  turning  angle,  whereas  the  maximum  turbulent 
boundary -layer  turning  angle  is  slightly  less  than  a  factor  of  three 
greater  than  the  corresponding  inviscid  value.  These  findings  again 
confirm  the  statement  by  Vaglio-Laurin  (Ref.  11)  quoted  earlier  that 
"due  to  the  larger  shearing  stress,  smaller  three-dimensional  effects 
can  be  expected  for  turbulent  layers  as  compared  with  laminar  layers 
subject  to  the  same  boundary  conditions"  for  the  case  of  sharp  cones  at 
incidence  in  a  hypersonic  flow  under  cold  wall  conditions. 

A  knowledge  of  wall  temperature  effects  on  the  three-dimensional 
hypersonic  turbulent  boundary  layer  is  of  interest  for  the  extrapolation 
of  results  from  ground  test  (wind  tunnel)  measurements  to  the  flight 
condition.  Adams  (Ref.  54)  has  recently  examined  the  three-dimensional 
laminar  boundary  layer  on  sharp  cones  at  incidence  under  hypersonic 
conditions  and  found  that  the  surface  upwash  angle,  ws,  may  be  substan¬ 
tially  increased  between  cold  wall  (Tw/T0  -*0)  conditions  representative 
of  hypersonic  atmospheric  flight  and  hot  wall  (TW/TQ  «  0.  9)  conditions 
representative  of  force  testing  in  hypersonic  wind  tunnels.  Figure  22 
presents  the  same  type  of  information  for  the  present  three-dimensional 
hypersonic  turbulent  boundary  layer  on  the  7.  2-deg  half-angle  sharp 
cone  at  4.  0-deg  angle  of  attack.  Note  from  Fig.  22  that  maximum  turn¬ 
ing  in  the  present  three-dimensional  turbulent  boundary  layer  occurs  at 
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0  ss  110  deg  irrespective  of  wall  temperature  ratio  in  the  range  con¬ 
sidered.  A  value  of  TW/TQ  =  0.  8  is  representative  of  force  tests  in  the 
AEDC-VKF  Tunnel  B,  whereas  TW/TQ  =0.4  is  representative  of  heat- 
transfer  tests  via  the  thin-skin  technique;  a  value  of  Tw/T0  =  0.  025  is 
representative  of  hypersonic  atmospheric  flight  under  the  same  Mach 
number— Reynolds  number  combination  as  the  wind  tunnel.  Figure  22 
shows  that  the  maximum  calculated  surface  upwash  angle,  us,  is 
approximately  8  deg  for  the  cold  wall  flight  condition.  This  is  a  factor 
of  two  greater  than  the  maximum  calculated  inviscid  turning  angle  of 
we  %  4  deg  as  determined  from  the  Jones  analysis  (Refs.  39  and  43). 

For  the  hot  wall,  force  test  condition  the  maximum  calculated  surface 
upwash  angle  is  approximately  13  to  14  deg,  a  factor  of  three  greater 
than  the  maximum  calculated  inviscid  turning  angle  of  ue  «  4  deg  as 
determined  from  the  Jones  analysis  (Refs.  39  and  43)  and  a  factor  of 
almost  two  greater  than  the  maximum  cold  wall  flight  value.  This 
means  that  for  the  present  sharp  cone,  angle  of  attack,  and  flow  condi¬ 
tions  the  maximum  surface  upwash  angle,  us,  of  the  three-dimensional 
turbulent  boundary  layer  may  be  increased  approximately  70  to  80  per¬ 
cent  between  cold  wall  flight  conditions  and  hot  wall  hypersonic  wind 
tunnel  force  test  conditions.  Although  this  difference  in  surface  turning 
angle  between  hot  and  cold  wall  conditions  for  the  three-dimensional 
turbulent  boundary  layer  is  not  as  large  as  the  differences  found  by 
Adams  (Ref.  54)  for  the  three-dimensional  laminar  boundary  layer, 
these  findings  indicate  to  the  vehicle  designer  that  proper  interpretation 
of  hot  wall  ground  test  measurements  on  slender  bodies  at  incidence 
under  hypersonic  conditions  with  either  a  laminar  or  turbulent  boundary 
layer  is  necessary  relative  to  cold  wall  flight  conditions  for  aerodynamic 
parameters  where  the  boundary -layer  flow  direction  is  important. 

Details  of  the  various  three-dimensional  turbulent  boundary -layer 
profiles  ait  the  $  =  90.0-deg  circumferential  location  and  x/L  =  0.  8  body 
station  for  the  present  sharp  cone  at  4.  0-deg  an^Le  of  attack  are  pre¬ 
sented  in  Fig.  23  relative  to  wall  temperature  effects.  Note  from 
Figs.  23a  and  b  that  in  the  near-wall  region  (y/L  <  10"^)  of  the  three- 
dimensional  turbulent  boundary  layer,  the  colder  the  wall,  the  larger 
the  numerical  value  for  both  the  streamwise  and  crossflow  velocity  at 
a  fixed  y/L  location.  Furthermore,  Fig.  23  c  shows  that  increasing  the 
wall  cooling  will  increase  the  gas  density  near  the  surface  (recall  that 
pa  1/T,  since  the  static  pressure  is  constant  across  the  boundary  layer). 
Hence,  the  effects  of  compressibility  on  the  three-dimensional  turbulent 
boundary  layer  are  such  as  to  result  in  increased  momentum  in  the 
lower  portions  of  the  boundary  layer  as  the  wall  cooling  is  increased. 
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Such  is  reflected  in  greater  turning  for  the  hot  wall  boundary  layer  rela¬ 
tive  to  the  cold  wall  case,  as  shown  in  Fig.  23d.  Further  observe, 
from  Fig.  23b,  that  in  the  outer  regions  of  the  three-dimensional  turbu¬ 
lent  boundary  layer,  the  hotter  the  wall,  the  larger  the  amount  of  cross - 
flow  velocity  profile  overshoot  relative  to  the  inviscid  crossflow  velocity 
value.  This  again  is  reflected  in  Fig.  23d  where  the  -hot  wall  boundary- 
layer  turning  is  greater  than  the  cold  wall  case.  As  shown  in  Fig.  23c, 
the  hot  wall  {Tw/T0  =  0.  8)  condition  representative  of  hypersonic  wind 
tunnel  force  testing  results  in  essentially  an  adiabatic  wall  condition, 
i.  e. ,  (3T/8y)  «  0.  Little  effect  is  seen  in  Fig.  23e  on  the  maxi- 

y 

mum  value  of  the  eddy- viscosity  ratio  with  respect  to  wall  temperature. 
However,  keep  in  mind  that  the  maximum  eddy  viscosity,  e  ,  must  in¬ 
crease  with  increasing  wall  temperature,  since  the  laminar  viscosity, 
ju,  increases  as  the  wall  temperature  is  increased,  with  the  maximum 
e/ju  ratio  remaining  constant.  All  of  the  above  findings  indicate  that,  in 
general,  the  hotter  the  wall,  the  greater  the  turning  effect  on  the  hyper¬ 
sonic  three-dimensional  turbulent  boundary  layer  attributable  to  cross - 
flow.  The  prior  work  by  Adams  {Ref.  54)  found  exactly  the  same  trend 
for  the  hypersonic  three-dimensional  laminar  boundary  layer. 


SECTION  V 

CONCLUDING  SUMMARY 


The  results  of  the  present  investigation  indicate  that  numerical 
calculation  of  the  three-dimensional  compressible  turbulent  boundary 
layer  on  a  sharp  cone  at  incidence  in  a  supersonic  or  hypersonic  stream 
is  indeed  feasible  and  reasonable,  based  on  comparisons  with  experi¬ 
mental  measurements.  The  assumption  of  a  locally  similar  turbulent 
boundary -layer  analysis  neglecting  "upstream  history"  on  the  sharp 
cone  appears  to  result  in  acceptable  calculations  for  the  mean  flow  pro¬ 
files,  including  crossflow,  when  used  in  conjunction  with  the  three- 
dimensional  invariant  turbulence  scalar  eddy-viscosity  model.  The 
main  advantage  of  the  present  approach  over  the  previously  used  integral 
techniques  is  that  no  assumptions  as  to  the  form  or  character  of  the  pro¬ 
files  are  required;  the  present  profiles  are  the  result  of  numerical 
integration  of  the  complete  governing  three-dimensional  turbulent 
boundary -layer  equations,  with  the  only  assumptions  being  the  three- 
dimensional  eddy -viscosity  formalism,  the  choice  of  mixing-length 
model,  and  the  locally  similar  method  of  solution.  The  degree  of 
success  experienced  in  the  present  investigation  indicates  that  the  invar¬ 
iant  turbulence  scalar  eddy-viscosity  approach  should  be  applicable  to 
numerical  calculation  of  general  three-dimensional  turbulent  boundary- 
layer  flows.  However,  much  more  work,  especially  carefully  controlled 
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experimental  investigations  of  the  three-dimensional  compressible 
turbulent  boundary-layer  structure,  remains  to  be  done  before  the 
question  of  variance  or  invariance  can  be  completely  resolved. 
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7.2-deg  Half-Angle  Sharp  Cone 
M,j,  *  8.0,  Re^/ft  -  3. 78x  10& 
Tw/T0  •  0.40,  L  •  3.66  ft 


of  Martellucci  and  Neff  (Ref.  47) 

Fig.  2  Spatial  Distribution  of  End  of  Transition  with  Respect  to  Angle  of  Attack 
under  Hypersonic  Conditions 
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12. 50-deg  Half-Angle  Sharp  Cone  at  a  ■  25. 0  deg 
Mn,  ■  1. 80,  Re00>  L  ■  2. 56  x  107,  Tw/T0  =  1. 0 
L  =  41. 58  In. ,  Turbulent  Boundary  Layer 


Fig.  3  Experimentally  Determined  Surface  Flow  Angles  at  Various  Lengthwise 
Stations  from  Rainbird  {Ref.  45) 
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12. 50-deg  Half-Angle  Sharp  Cone  at  a  ■  15. 780  deg 


Fig.  4  Surface  Pressure  Distribution  under  Supersonic  Conditions 
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12.50-deg  Half-Angle  Sharp  Cone  at  a  -  15.780  deg 
Ma,  “1.80,  y  - 1.40 


I  nviscid  Sharp  Cone  at  Angle  of  Attack 
Following  Jones  (Refs.  39,  43) 


0,  deg 

Fig.  5  Calculated  Inviscid  Flow  Parameters  on  Cone  Surface  under  Supersonic  Conditions 
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u/Ue 


a.  Streamwise  Velocity  Ratio 

Fig.  6  Windward  Ray  Turbulent  Boundary- Layer  Profiles  under  Supersonic 
Conditions 


46 


b.  Static  Temperature  Ratio 
Fig.  6  Concluded 
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^Ue 

a.  Streamwise  Velocity  Ratio 

Fig.  7  Comparison  of  Windward  Ray  Turbulent  Boundary- Layer  Profiles  with 
Previous  Analysis  by  Adams  (Ref.  48) 
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b.  Static  Temperature  Ratio 
Fig.  7  Concluded 
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y.  in. 


12. 50-deg  Half-Angle  Sharp  Cone  at  a  =  15. 780  deg 
Mm  - 1 80,  Re^  l  =  2. 56  x  107.  Tw/T0  -  L  0 
x/L  =  0. 85,  L  -  41. 58  in. ,  0  =  45. 0  deg 


- Three-Dimensional  Turbulent  Boundary-Layer  Theory 

g  o  A  Experimental  Data  from  Rainbird  (Ref.  45) 


a.  0  =  45.0  deg 

Fig.  8  Three-Dimensional  Turbulent  Boundary-Layer  Profiles  under 
Supersonic  Conditions 
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12.50-deg  Half-Angle  Sharp  Cone  at  a  =  15. 780  deg 
Mp  =  1. 80,  Reroi  L  =  2. 56  x  107,  yT0  - 1. 0 

x/L  =  0.85,  L  3  41.58  in.,  0  =  90.0  deg 

- Three-Dimensional  Turbulent  Boundary-Layer  Theory 

a  o  A  Experimental  Data  from  Rainbird  (Ref.  45) 


b.  <p  =  90.0  deg 
Fig.  8  Continued 
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0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8  2.0 


c.  0  =  135.0  deg 
Fig.  8  Concluded 


12.50-deg  Half-Angle  Sharp  Cone 
a  -  15. 780  deg,  Mm » 1. 80 
Rem  ■  «  2.56x10';  W  L0 
rfL-ft*  L-4L581P. 


Fig.  9  Streamline  Directions  within  the  Three-Dimensional  Turbulent  Boundary  Layer 
under  Supersonic  Conditions 
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12.50-deg  Half-Angle  Sharp  Cone  at  a  ■  15. 780  deg 
- 1. 80,  Rem  ,  -  2.56  x  10?,  TJT0  *  1.0 

x/L-0.85,  L-  41.58  in. 

- Three-Dimensional  Laminar  Boundary-Layer  Theory 

...  Three-Dimensional  Turbulent  Boundary-Layer  Theory 
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Fig.  1 1  Eddy-Viscosity  Distributions  at  Various  Angular  Locations  under 
Supersonic  Conditions 
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12. 50-deg  Half-Angle  Sharp  Cone  at  a  =  15. 780  deg 
Mqo  =  1. 80,  RCqo  l  ■  2. 56  x  10?  Tw/T0  =  1. 0 
x/L  =  a 85,  L  =  41. 58  in.,  0  -90.0 deg 

- Three-Dimensional  Turbulent  Boundary-Layer  Theory 

- Three-Dimensional  Laminar  Boundary-Layer  Theory 


a.  u  -  Velocity  Ratio 

Fig.  12  Comparison  of  Laminar  and  Turbulent  Three-Dimensional 
Boundary- Layer  Profiles  at  Angular  Location  <f>  =  90.0  deg 
under  Supersonic  Conditions 
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12. 50-deg  Half-Angle  Sharp  Cone  at  a  =  15. 780  deg 
Mqq  - 1. 80,  Rem  L  -  2. 56  x  107,  Tw/T0  -  1. 0 
x/L  ■  0. 85,  L  -  4L  58  in. ,  0  ■  90. 0  deg 

- -  Three-Dimensional  Turbulent  Boundary-Layer  Theory 

— - Three-Dimensional  Laminar  Boundary-Layer  Theory 


1.0  1.1  1.2  L3  L4  1.5  1.6 

T'Te 

c.  Static  Temperature  Ratio 
Fig.  12  Concluded 
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«e 

■  I  .I  i  Ug)  i^=^=t=  i  i  i  1  i 

0  20  40  60  80  100  120  140  160  180 

0,  deg 

a.  a  =  1 .0  deg 

Fig.  13  Calculated  Inviscid  Flow  Parameters  on  Cone  Surface  under 
Hypersonic  Conditions 
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b.  a  =  2.0  deg 
Fig.  13  Continued 
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c.  a  =  4.0  deg 
Fig.  13  Concluded 
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7. 2-deg  Half-Angle  Sharp  Cone  at  a  -  4. 0  deg 

Mq,  -  7. 97,  Re^/ft  -  1. 32  x  106,  Tw/T0  •  0.42,  L  -  3.66  ft 


0  20  40  60  80  100  120  140  160 


0.  deg 

Fig.  14  Laminar  Heat-Transfer  Distribution  on  a  Sharp  Cone  at  Incidence  under  Hypersonic 
Conditions 
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x/L 


a.  a  =  1.0  deg 

Fig.  15  Turbulent  Heat-Transfer  Distribution  on  a  Sharp  Cone  at 
Incidence  under  Hypersonic  Conditions 
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7.2-deg  Half -Angle  Sharp  Cone  at  a  ■  2.0  deg 
Mq)  ■  8. 0,  Reoo/ft  -3.”78x  ltf 


b.  a  =  2.0  deg 
Fig.  15  Continued 
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- )( - Three-Dimensional  Turbulent  Boundary -Layer  Theory 

o  A  b  o  □  Experimental  Data  from  Martellucci  and  Neff  {Ref.  47) 
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Fig.  16  Circumferential  Turbulent  Heat-Transfer  Distribution  on  a  Sharp  Cone 
at  Incidence  under  Hypersonic  Conditions 


c  c 


ili; 


Inner  Region  (y/6t  <  0.2) 
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Fig.  17  Windward  Ray  Turbulent  Boundary- Layer  Profiles  Under  Cold  Wall  Hypersonic 
Conditions 
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Fig.  18  Windward  Ray  Turbulent  Boundary-Layer  Parameters  under  Cold  Wall  Hypersonic 
Conditions 
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b.  Crossflow  Velocity  Ratio 
Fig.  20  Continued 
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c.  Static  Temperature  Ratio 
Fig.  20  Continued 
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d.  Streamline  Directions  in  Boundary  Layer 
Fig.  20  Continued 
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e.  Eddy-Viscosity  Ratio 
Fig.  20  Concluded 
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7.2-deg  Half-Angle  Sharp  Cone  at  a  -4.0  deg 
Mm-8.a  Rem/ft-3.78xlOP 

T^-0.40,  L- 3.66  ft,  x/L-0.80 

— Three-Dimensional  Laminar  Boundary-Layer  Theory 
—  Three-Dimensional  Turbulent  Boundary-Layer  Theory 


Fig.  21  Surface  and  External  Flow  Directions  under  Cold  Wall  Hypersonic  Conditions 
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7. 2-deg  Half-Angle  Sharp  Cone  at  a  -  4. 0  deg 
Mq,  •  8. 0;  Rejft  ■  3.78  x  10P.  x/L  -  0. 80,  L  -  3.66  ft 


Fig.  22  Effects  of  Wall  Temperature  on  Surface  Flow  Direction  under  Hypersonic  Conditions 
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a.  Stream  wise  Velocity  Ratio 

Fig.  23  Effects  of  Wall  Temperature  on  Three-Dimensional  Turbulent 
Boundary- Layer  Profiles  under  Hypersonic  Conditions 
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b.  Crossflow  Velocity  Ratio 
Fig.  23  Continued 


■  4.0  deg 

o.sa 
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0123  456  7  89  10  11  12  13 

u.  deg 

d.  Streamline  Directions  in  Boundary  Layer 
Fig.  23  Continued 
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TABLE  I 

COMPARISON  OF  EXACT  RESULTS  WITH  APPROXIMATE  SCALING  TECHNIQUE  FOR  STANTON  NUMBER 
DISTRIBUTIONS  ON  A  SHARP  CONE  AT  INCIDENCE  UNDER  COLD  WALL  HYPERSONIC  CONDITIONS 


oo 

to 


st. 

x  103 

x/L 

4  = 

0  deg 

4  = 

45  deg 

4  = 

90  deg 

4  = 

135  deg 

4  = 

180  deg 

Exact 

Eq.  (49) 

Exact 

Eq.  (49) 

Exact 

Eq.  (49) 

Exact 

Eq.  (49) 

Exact 

Eq.  (49) 

0.50 

1.  959 

1.  981 

1.  710 

1.  726 

1.204 

0.824 

0.833 

ng 

0.  655 

0.  65 

1.867 

1.879 

1.  629 

1.  638 

1.  149 

0.  786 

0.  790 

RHH 

0.  621 

0.  80 

1.803 

1.803 

1.  571 

1.  571 

1.  105 

1.  105 

0.  758 

0.  758 

0.596 

0.  596 

0.  95 

1.  748 

1.  742 

1.  526 

1.  518 

1.074 

1.068 

0.  735 

0.  732 

0.  580 

0.  576 

7.  2 -deg  Half -Angle  Sharp  Cone  at  a  =  4.  0  deg 
M,,,  =  8.  0,  Re„/ft  =  3.  78  x  106 
Tw/T0  =  0.40,  L  =  3.  66  ft 

Note:  In  the  above  table,  "exact"  denotes  the  present  three-dimensional 
turbulent  boundary-layer  theory. 
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